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The  state  and  control  derivatives  of  the  identified  model  agree  well 
witii  the  wind  tunnel  model.  Most  of  the  identified  dynamic  derivatives 
compare  well  with  the  theoretically  predicted  value  at  low  alpha  and  some 
compare  well  at  high  alpha.  Qne  identified  dynamic  derivative  which 
doesn't  compare  well  at  hign^^ha  is  .  At  high  alpha,  is 

predicted  to  have  values  above  ^0;  the  ic^ntified  model  gives  ^lues 
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LIST  OP  SYMBOLS 


Variable 


Description 


Units 


Longitudinal  acceleration  measured  along  (m/sec^) 

aircraft  x-body  axis  at  L  i  ■  nose 

Lateral  acceleration  measured  along  aircraft  (m/sec^) 

y-body  axis  at  i  *  eg.  nose*  tail 

Vertical  acceleration  measured  along  aircraft  (m/sec^) 

z-body  axis  at  i.  i  •  eg,  nose*  tail*  R*  L 

Reference  wingspan  (m) 

Roll  moment  coefficient  about  x*body  axis 
Pitch  moment  coefficient  about  y>body  axis 
Yaw  moment  coefficient  about  z>body  axis 
Axial  force  coefRcient  along  x-body  axis 
Side  force  coefficient  along  ybody 

Normal  force  coefficient  along  z-body  axis 

6  Cy/  ip  t  Other  comMnations  occur  where  y  is 
changed  to  x*  z*  i  *  m*  n  and  p  is  changed  to 
another  state,  to  a  control,  or  to  ff  or  0 . 

Normal  force  coefficient  along  z>body  axis 

Mean  aezxxfynaxnic  chord  (m) 

Acceleration  due  to  gravity  «  9.806  m/sec^  (m/sec^) 

Angular  momentum  at  both  T-2C  engines  (Newton-m-sec) 


Moment  of  inertia  about  x  body  axis 
Moment  of  inertia  about  y  body  axis 
Moment  of  inertia  about  z  body  axis 
Cross  product  of  inertial  about  x,  z  body  axes 


(kg  m^) 
(kgm^ 
(kg  m*) 
(kg  m^) 


Variable 

Description 

Units 

^x 

Distance  along  x-axis  between  aerodynamic 
center  and  center  of  gravity  (positive  for 
a.  c.  forward  of  c.  g. ) 

(m) 

Distance  along  z-axis  between  aerodynamic 
center  and  center  of  gravity  (positive  for 
a.  c.  below  c.  g. ) 

(m) 

"z 

e 

Distance  along  z-axis  between  thrust  axis 
and  center  of  gravity  (positive  for  thrust 
axis  below  c.  g.  ) 

(m) 

m 

Aircraft  mass 

(kg) 

P 

Aircraft  x-body  axis  roll  rate 

(rad/ sec) 

q 

Aircraft  y-body  axis  pitch  rate 

(rad/ sec) 

r 

Aircraft  z-body  axis  yaw  rate 

(rad/ sec) 

S 

Reference  wing  area 

(m^) 

T 

X 

Thrust  component  silong  x-body  axis  direction 

(Newtons) 

T 

z 

Thrust  component  along  z-body  axis  direction 

(Newtons) 

u 

Speed  along  aircraft  x-body  axis  with  respect 
to  a  nonmoving  free  stresim  air  mass 

(m/sec) 

V 

Aircraft  total  airspeed 

(m/ sec) 

V 

Speed  along  aircraft  y-body  axis  with  respect 
to  a  nonmoving  free  stream  ear  mass 

(m/sec) 

w 

Speed  along  aircraft  z-body  axis  with  respect 
to  a  nonmoving  free  stream  air  mass 

(m/sec) 

XV 


Description 


Units 


Aircraft  c.  g.  ;position  along  the  longitudinal  (nx) 

axis  of  earth  fixed  axes  oriented  with  respect 
to  aircraft  body  axes  by  Euler  angles.  Positive 
when  c.  g.  is  forward  of  axis  system  origin. 

Distance  between  c.g.  and  sensor  i  measured  (m) 
along  x-body  axis,  .  V.  Positive  for  sensor 

forward  of  c.  g. 

Distance  between  c.g.  and  longitudinal  accelero-  (m) 
meter  at  i  measured  along  x>body  axis, 
i  «  NOSE.  Positive  for  accelerometer  forward 
of  c.  g. 

Distance  between  c.  g.  and  lateral  accelerometer  (m) 
at  i  measured  along  x-body  axis,  i  >  c.g. ,  nose, 
tall.  Positive  for  accelerometer  forward  of  c.  g. 

Distance  between  c,  g.  and  vertical  accelerometer  (m) 
at  1  measured  along  x-bo47  axis,  i  >  c.g.  nose 
tail.  Positive  for  accelerometer  forward  of 
c.g. 

Aircraft  c.  g.  position  along  the  lateral  axis  of  (m) 
earth  fixed  axes  oriented  with  respect  to  air¬ 
craft  body  axes  oriented  with  respect  to  aircraft 
body  axes  by  Euler  angles.  Positive  when  c.  g. 
is  right  of  axis  system  origin. 

Distance  between  c.  g.  and  sensor  i.  Positive  for  (m) 
sensor  right  of  c.g. 

Distance  between  c.  g.  and  longitudinal  accelero-  (m) 
meter  at  i  measured  along  y-body  axis,  i  »  NOSE, 
Positive  for  accelerometer  right  of  c.  g. 
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Variable 


Description 


Units 


Distance  between  c.g.  and  lateral  accelero* 
meter  at  i  measured  along  y-body  axis. 

1  «  c.  g. «  NOSE.  TAIL.  PositiTe  for 
accelerometer  rif^t  of  c.  g. 

(m) 

Distance  between  c.g.  and  vertical  accelero- 
meter  at  i  measured  along  ybody  axis, 
i  -  CG.  NOSE.  TAIL.  R.L.  Positive 
for  accelerometer  right  of  c.  g. 

(m) 

z 

Aircraft  c.g.  position  aloi^  the  vertical  axis 
of  earth  fixed  axes  oriented  with  respect  to 
aircraft  body  axes  by  Euler  angles.  Positive 
when  c.g.  is  below  axis  system  origin. 

(m) 

Distance  between  c.g.  and  sensor  i  measured 
along  z-body  axis,  i  *  a.  d.  V.  Positive  for 
sensor  below  c.g. 

(m) 

“z 

*l 

Distance  between  c.  g.  and  longitudinal  accelero* 
meter  at  i  measured  along  z-body  axis, 
i  ■  NOSE.  Positive  for  accelerometer  below 

C« 

(m) 

z 

Distance  between  c.  g.  and  lateral  accelero¬ 
meter  at  i  measured  along  z-body  axis, 
i  ■  CG.  NOSE.  TAIL.  Positive  for  accelero¬ 
meter  below  c.g. 

(m) 

*l 

Distance  between  c.  g.  and  vertical  accelero¬ 
meter  at  i  measured  along  z-body  axis. 

1  ■  c.g..  NOSE.  TAIL,  R.L.  Positive  for 
accelerometer  below  c.  g. 

(m) 
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Variable 


Description 


Dmts 


a 

iiagle  of  attack 

(rad) 

• 

a 

Time  derivative  of  a 

(rad/ see) 

Sidesl4>  angle 

(rad) 

• 

P 

Time  derivative  of 

Crad/sec) 

Aileron  de£leetioa(  +  5  -•  left  wing  down  roll) 

(rad) 

ft 

e 

Elevator  deflection (^f  5^  ^  nose  down  pitch) 

(rad) 

Rudder  deflection (  +  6^  -  nose  left  yaw) 

(rad) 

e 

Euler  pitch  angle 

(rad) 

p 

Air  density 

(kg/m®) 

P 

Euler  bank  axigle 

(rad) 

• 

LIST  OF  SUBSCRIPTS 

o 

Denotes  reference  value 

c,g. 

Center  of  gravity 

L 

Left  wlngtip 

m 

Designates  a  measured  quantity  when  used  as  a 
subscript 

R 

Right  wingtip 

a 

Stability  axes 

NOSE 

Nose 

TAIL 
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1.  INTRODUCTION 

1.1  BACKGROUND 

Our  previous  effort  which  is  documented  as  Volume  n  (Ref. 

[1])  has  demonstrated  that  the  Estimation-Before-Modeling  (EBM) 
system  identification  method  extracts  with  precision  the  nonlinear 
aerodynamics  of  the  high  angle  of  attack/ sideslip  domain  for  aircraft. 

This  demonstration  was  carried  out  using  synthetically  generated  data 
from  a  T-2C  wind  tunnel  model  and  a  set  of  theoretically  predicted 
d3mamic  derivatives,  provided  by  die  Naval  Air  Development  Center 
(NADC).  The  data  of  that  simulation  study  were  generated  by  exciting 
the  synthetic  model  with  the  controls  and  the  initial  states  contained  in 
actual  T-2C  flight  data  provided  by  NADC.  The  T-2C  light  jet  trainer 
aircraft  was  instrumented  and  operated  by  NADC  for  the  express  purpose 
of  providing  high  angle  of  attack/ sideslip  flight  data  for  system  identifica¬ 
tion  studies.  The  measurement  equations  and  the  noise  levels  of  the  col¬ 
lected  data  were  used  in  generating  the  synthetic  measurements  for  the 
simulation  study.  As  a  result,  the  EBM  methodology  was  tested  under 
controlled  but  realistic  simulated  conditions.  The  identified  model 
closely  follows  the  nonlinearities  of  the  dynamic  derivatives  and  is  sm 
excellent  match  to  the  control  and  state  derivatives  of  the  wind  tuzmel 
model.  The  predicted  responses  of  the  identified  model  are  good  matches 
to  the  corresponding  responses  of  the  wind  tunnel  model.  The  work  docu¬ 
mented  in  Volume  II  validates,  therefore,  the  excellent  performance 
capability  of  the  EBM  system  identification  method.  This  present  volume 
contains  the  application  of  this  method  to  real  data. 

1. 2  SCOPE  OF  FLIGHT  DATA  ANALYSIS  STUDY 

The  EBM  system  identification  method  is  used  to  process  eighteen 
T-2C  flight  maneuvers  containing  over  600  seconds  of  data  sampled  at 
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20  Hertz.  The  instruxaentation  of  the  T-2C  provided  24  measured 

quantities  including  three  controls.  The  first  step  of  the  EBM 

method  provides  optimal  estimates  of  the  adrcraft  states,  of  the 

aerod3mamic  forces  and  moments  and  of  the  biases  and  scale  factors 

in  the  measurements.  The  second  step  of  the  EBM  method 

identifies  detailed  models  of  the  lateral  aerodynamic  coefficients  C 

and  and  of  the  longitudinal  aerodynamic  coe^cient  C^. 

Time  did  not  permit  the  identification  of  a  detail  model  of  and 

C  :  but  partial  models  are  presented.  The  identified  models  are 
z 

compared  to  the  T-2C  wind  tunnel  model  and  to  the  theoretical 
predictions  of  Ref.  [2]. 


f 


1.3  ESTIMATION  BEFORE  MODELING  (EBM)  METHOD 

The  Estimation  Before  Modeling  (EBM)  system  identification 
method  is  a  two -step  approach  to  the  modeling  of  aircraft  aerodynamics. 
The  first  step  uses  Gauss -Markov  processes  to  model  the  accelerations 
due  to  aerodynamic  forces  and  moments  and  it  employs  an  extended 
Kalman-Bucy  filter/Bryson-Frazier  smoother  to  generate  estimates  of 
states,  aerodynamic  forces  3md  moments  time  histories,  and  measure¬ 
ment  biases  and  scale  factors.  A  description  of  the  filter /smoother  is 
presented  in  Appendix  B  of  Volume  n  (Ref.  [1]). 

In  the  second  step.  Stepwise  Multiple  Linear  Regression  (SMLR) 
is  used  together  with  subspace  modeling  to  identify  the  dependency  of  an 
aerodynamic  coefficient  on  aircraft  states  and  controls.  A  detailed 

account  of  step  two  is  presented  in  Section  1. 3  and  Appendix  C  of 
Volume  II  (Ref.  [1]. 


2.  SCENAEIO  OF  THE  T-2C  SYSTEM  IDENTIFICATION  PROBLEM 


2. 1  FLIGHT  TEST  CONDITIONS  AND  DESCRIPTIONS  OF 
MANEUVERS 

The  T-2C  flight  test  data  were  gathered  with  gear  up  and 
speed  brake  retracted  under  the  following  conditions: 

True  Airspeed  -  75  to  150  meters/sec 

Mach  No.  •  0.25  to  0.4 

Altitude  -  6000  to  7600  meters 

The  throttles  were  fixed  during  maneuvers.  The  Naval  Air  Develop¬ 
ment  Center  provided  fovir  sets  of  flight  data  consisting  of  eighteen 
maneuvers,  amoimting  to  over  600  seconds  of  20  Hertz  data.  The 
type  of  input  used  to  generate  each  maneuver  is  described  in  the  last 
column  of  Table  2.1-1.  The  first  coluixm  gives  the  date  of  the  flight 
set  and  the  second  column  denotes  the  flight  number  and  the  maneuver 
number  that  goes  with  the  input  used. 

The  time  histories  of  the  inputs  are  presented  in  Appendix  B 
for  all  of  the  maneuvers.  The  maneuvers  FlMl,  F4M9  and  F6M1  have 
360°  rolls  and  the  maneuver  F4M10  has  two  360®  rolls.  The  angle  of 
attack,  a,  ranges  from  -4®  to  40°  euid  the  sideslip  angle,  P,  covers 
tile  range  from  -14*  to  26*.  Most  of  the  maneuvers  have  a  trim  a  be¬ 
tween  10. 5°  and  11.5*  and  a  trim  V  (i.  e. ,  true  airspeed)  between  90.0 
and  95.0  meters/sec.  The  maasuvers  F4M7  and  F4M8  have  a  trim  °- 
between  5.5®  and  6.5®  and  a  trim  V  between  115  and  120  meters/sec. 


TABLE  2. 1-1  Description  of  the  Actual  T-2C  Data  Processed 
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2.  2  T-2C  INSTRUMENTATION  -  MEASUREMENT  EQUATIONS 

Table  2.  2-1  lists  the  twenty -four  variables  that  are  sensed 
by  the  T-2C  instrumentation  package  and  provides  measurement  range 
and  location  with  respect  to  the  normal  flight  test  center  of  gravity  for 
each  sensor.  This  package  as  well  as  the  data  acquisition  techniques 
are  essentially  the  same  as  those  described  in  detail  in  Ref.  [3]  for 
the  NADC  YT-2B  airplane  and  will  not  be  discussed  here.  Three  of 

the  twenty-four  variables  are  control  variables  and  three  others  are 

•  •  • 

the  angular  acceleration  variables  p,  q  and  r.  The  angular  accelera¬ 
tion  measurements  were  not  of  good  quality  and.  therefore,  were  not 
used  in  our  study. 

The  measurement  equations  relating  the  sensed  angle  of  attack, 
sideslip  angle  and  total  airspeed  variables  to  aircraft  body-axes  states 
are: 

-1  /w-qS^tpy, 

°  \  u  +  q  z  “  r  y 

a  ct 

,  I  V  -  p  z  -  +  r  Xo 
m  'u  +  qz^-ry^/ 

=  [(u  +  q  iy  -  r  y^)  +(v  +  r  Xy  -  p  Zy)  +  (w  +  p  yy  -  q  Xy)  V 

The*  measurement  equations  relating  the  sensed  roll,  pitch  and 
yaw  rates,  attitude,  altitude,  angular  accelerations  smd  control  variables 
to  aircraft  states  are  as  follows: 
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The  nine  measurement  equations  relating  the  sensed  linear 
acceleration  variables  to  aircraft  body-axes  states  are  as  follows: 


B  -  (r^  +  q^)  X  +  (p  q  -  r)  y„  +  (p  r  +  q)  z 

1  Xf  X.  X. 


•  •  2  2 

B,  +  (p  q  +  r)  X  -  (p  +  r  )  y  +  (q  r  -  p)  z 


•  _  -  _  2  2  — 

B_  +  (p  1'  -  q)  X  +  (q  r  +  p)  y  -  (p  +  q  )  z 

^  \  \  \ 


where  i  -  nose,  j  =  c.g.,  NOSE  and  TAIL,  and  k  =  c.g. ,  NOSE, 

TAIL,  R  and  L.  The  values  of  the  x's,  y's  and  z's  are  given 

in  Table  2.2-1.  For  example,  x  =4.62  meters, 

*NOSE 

y  =  -  .  13  meters  and  z  =  0. 18  meters.  The  states 

*NOSE  *NOSE 

B^,  B2  and  B^  are  the  fore-aft,  the  lateral,  and  the  normal  accel 
erations  at  c.g.,  respectively. 


Table  2.  2-1  Location  Values  of  instruments 
in  ^e  T-2C  Instrumentation  Package; 

Positive  z  means  Instrument  is  Foward  of 

c.g. .  Positive  y  •  right  of  c.g..  Positive 

2  •  below  c.g. 

Measured  Quantity 

Location  in  Meters 

Measurement 

No. 

£ 

f 

z 

Range 

1. 

i^ngle  oi  Attack  (a  ) 
ni 

7.04 

0 

0 

-15®  to  45“ 

2. 

Sideslip  Angle  (S  ) 
m 

7.13 

0 

0 

±  30® 

3. 

Airspeed  (V  ) 
m 

4.62 

0.13 

0.18 

0-500  kt;  ±  50  kt 

4. 

Boll  Rate  (p  ) 
m 

4.04 

0.13 

0.05 

+  ISO  deg/sec 

5. 

Pitch  Rate  (q^) 

0.12 

0 

0.58 

+  45  deg/sec 

6. 

Yaw  Rate  (r  ) 
od 

4.04 

-.13 

0.05 

150  deg/ sec. 

7. 

Pitch  Angle  (6  ) 
m 

4.04 

0.13 

0.05 

^  90  deg. 

8. 

Roll  Angle  (i  ) 
ni 

4.04 

0.13 

0.05 

+  180  deg. 

9. 

Altitude  (h  ) 
nx 

4.04 

0.13 

0.05 

0  -  10670  m 

10. 

Roll  acceleration  (p^^) 

0.12 

0 

0.58 

^  275  deg/sec. 

11. 

Pitch  acceleration  (q^ 

0.12 

0 

0.58 

+  450  deg/aec. 

12. 

Yaw  acceleration  (r  ) 

01 

4.62 

0 

0 

+  275  deg/sec. ^ 

13. 

Fore-Aft  Accleration 

*oi  NOSE 

4.62 

-0.13 

0.18 

+  1.0  g 

14. 

15. 

Lateral  Acceleration  at  c.g. 

(ay  ) 

•^m  c.g. 

Lateral  Acceleration  at  NOSE 

0.12 

-0.02 

0.58 

+  1.0  g 

^  NOSE 

4.  62 

-0.13 

0.18 

+  l.Og 

16. 

Later^  Acceleration  at  Tail 
(a.  ) 

TAIL 

-3.89 

0.71 

-0.08 

+  l.Og 

17. 

Normal  Acceleration  at  c.g. 

m  c.g. 

0.12 

0.02 

0.58 

+  3.5g 

18. 

Normal  Acceleration  at  Nose 

*ni  NOSE 

4.04 

-0.13 

0.05 

+  5.0g 

19. 

Normal  Acceleration  at  TAIL 

*ni  TAIL 

-3.89 

0.71 

-0.08 

+  3.5g 

20. 

Normal  Acceleration  at  Right 
Wing  Tip  (a*^)  n 

-0.89 

5.08 

0.05 

+  3.5g 

i 


2-5 


Table  2.2-1 

(Continued) 

Location  in  Meters 

Measurement 

No. 

Measured  Quandtf 

S 

7 

r 

Range  , 

21 

Normal  Acceleration  at  Left 

-0.89 

-5.08 

0.05 

+  3.5  g 

Wing  Tip  <a,  ) 

*m  L 

22 

Aileron  (S^  ) 
m 

-0.97 

4.06 

0.05  . 

13  deg 

23 

Elevator  (t>  ) 

-5.01 

0.05 

-1.19 

-27.5  to  15  deg. 

®m 

24 

Rudder  (fi  ) 
r 
m 

-5.01 

0.05 

-0.72 

+  24  deg. 
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2.3  THE  T-2C  AIRCRAFT  CONSTANTS 

The  aircraft  constants  required  in  the  kinematical 
and  force  and  moment  equations  of  the  T-2C  are  given  in  this  section. 

V 

The  reference  wing  area,  the  mean  aerodynamic  chord  and 
the  reference  wingspan  are 

S  =  23.  7  m^ 

7  *  2.  26  m 

b  =  11.62  m 

The  distance  along  the  z-axis  between  the  aerod3niaznic  center 
and  the  center  of  gravity  and  the  distance  along  the  z-axis  between  the 
thrust  axis  and  the  center  of  gravity  are 

i  =  0  m 
z 

I  =  .49  m 
z 

e 

The  moments  of  inertia  about  the  body  axes  are 
I  =  12212  kg  m^ 

I  =  19811  kg  m^ 

y 

I  =  25781  kg  m^ 
z 

I  =  0  kg  m^ 
xz  ® 

The  total  angular  momentum  of  both  engines  is 

h  =  2448  Newtons-m-sec 
e 

X 

The  engines  do  not  covuiter  rotate.  They  rotate  clockwise  from  the 
rear  of  the  aircraft. 
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The  values  for  the  mass  m.  the  thrust  along  the  x  body- 

axis  and  the  moment  arm  e. ,  the  distance  along  the  x-axis 

betareen  the  aerodynamic  center  and  the  center  of  gravity)  are  given 

in  Table  2.3-1  The  thrust  T  along  the  z  body  axis  is  zero. 

z 
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Table  2. 3-1  Values  of  the  T-2C  Aircraft  Constants  that  Vary  from 
Maneuver  to  Maneuver 


Maneuver 

Number 

m 

Mass 

(kg) 

T 

Thrust 

(Newtons) 

i 

X 

Moment 

Arm 

(Meters) 

PlMl 

5275 

4130 

.042 

F2M1 

5330 

4092 

.065 

F2M2 

5198 

3968 

.072 

F2M3 

5039 

3884 

.054 

F3M1 

5384 

4119 

.055 

F3M2 

5039 

3823 

.046 

F3M3 

4944 

3835 

.068 

F4M1 

5420 

4107 

.041 

F4M2 

5371 

4154 

.044 

F4M3 

5325 

4012 

.021 

F4M4 

5266 

3990 

.021 

F4M5 

5216 

3991 

.  050 

F4M6 

5171 

3928 

;  .048 

F4M7 

5126 

3229 

.  203 

F4M8 

5080 

3226 

.177 

F4M9 

5035 

1 

3823 

.017 

F4M10 

1  4990 

:  3828 

.030 

F6M1 

1  5207 

] 

1  4486 

! 

.  006 

2.4  EQUATIONS  OF  MOTION 

The  following  kineznatical  relations  for  a  rigid  aircraft  are 
given  with  reference  to  body  axes  at  Ihe  center  of  gravity: 

# 

u®rv-qw-g  sin  'S  + 

• 

V  *  pw  -  r  u  +  g  cos  0  sin  0  +  Bg 
w®qu-pv  +  g  cos  6  cos  0  +  Bg 
;.pqC^j  +  qrC^2+,C^3  +  B^  + 

0  -  q  cos  0  -  r  sin  0 

0 

0  =  p  +  q  tan  0  sin  0  +  r  tan  0  cos  0 

0 

z  a  -  u  Sin  0  +  v  cos  0  sin  0  +  w  cos  0  cos  0 

where  B  ,  B  , B  are  the  accelerations  at  the  center  of  gravity 
12  6 

due  to  aerodynamic  and  thrust  forces  and  moments  acting  on  the  body 
of  the  aircraft. 

The  constant  coefficients  through  Cg^  are  given  by: 


"41 


I  I  -  I  , 
X  z  xz 


I  - 1 

Z  X 


"51 


C  * 
42 


I  (I  -  I  )  -  I 
z ' y  z  Xz 


I  I  -  I 
X  z  xz 


1 

xz 
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ESTIMATION  RESULTS 


The  extended  Kalman-Bucy  filter/Bryson-Frazier  smoother 
as  discussed  in  Ref.  [1]  was  used  to  process  the  ei^teen  maneuvers 
FlMl.  F2M1. ....  F6M1  to  obtain  the  estimation  results  of  this  study. 
Bias  states  and  scale  factor  states  were  added  to  the  filter /smoother 
to  handle  the  instrumentation  biases  and  scale  factors  that  are  inherent 
in  real  data. 

In  addition  to  using  die  ^ter/ smoother,  a  simple  technique 
was  used  in  determining  the  scale  factors  of  a  and  /3.  Let  I  be  the 
number  of  times  that  the  time  history  of  the  measured  a  crosses  its 
trim  value  anH  the  J  be  the  number  of  such  times  for  /9.  All  measured 
quantities  are  corrected  to  c.  g.  location.  The  v  and  w  equations 


-  r  u  +  g  cos  0  sin  0  +  a 


w  *qu-pv  +g  cos  0  cos  0  +  a  +  b, 

^  Bi  «  i 

c.g. 

with  u  =  y  -  (v^  +  w^)  are  integrated  using  pre-smoothed  values 
a  a 

of  V,  p,  q,  r,  0,  0,  a  and  a  ;  the  acceleration  biases  b^ 

^c.g.  ^c.g, 

and  b.  are  computed  over  the  intervals  t.]  and  [t.  t  ], 

respectively,  so  as  to  satisfy  the  state  constraints  a  (t.)  =  a(t.), 
i  »  1,  2, . . . ,  I  and  (L)  *  (L)  j  *  1,  2, ....  J  where  and 

are  the  integrated  values.  The  angle  of  attack  scale  factor  is  the 

& 

correlation  coefficient  of  the  two  time  histories  a  (t)  and  a  (t).  The 

a 

sideslip  angle  scale  factor  is  similarly  given. 


3. 1  COMPARISON  OF  ESTIMATED  VALUES  WITH  MEASURED 
VALUES 

The  filter /smoother  which  uses  the  equations  of  motion  of 

Section  2.4  provided  estimates  of  the  states  u,  v,  w,  p,  q,  r,  B,  0, 

altitude,  the  accelerations  . . . .  Bg  and  the  biases  and  scale 

factors.  These  estimated  values  were  used  in  the  measurement 

equations  of  Section  2. 2  to  obtain  the  curves  labeled  "Estimated 

Values"  that  are  plotted  in  Appendix  B.  Figures  B-1  through  B-18 

contain  the  results  for  each  of  the  eighteen  msmeuvers.  There  are 

four  pages  for  each  maneuver.  The  first  page  shows  the  plots  of 

a  and  /3.  The  second  presents  the  measured  values  of  the  controls 

6.6,6;  the  smoothed  values  of  the  controls  are  not  compared  in 
ear 

these  figures  with  the  measured  values.  The  third  page  consists  of 

the  variables  V.  p,  q,  r,  B,  0,  altitude  and  a  .  The  fourth 

z 

c.  g. 

page  covers  the  eight  accelerations  a  .  a  .  a  ,  a  , 

®R  ®L  ®NOSE  ®TAIL 

a  ,  a  .a  and  a  .  The  measured  values  of 

%OSE  ^TAIL  ^c.g.  *NOSE 

the  sensed  quantities  with  biases  and  scale  factors  removed  are  plotted 

using  a  solid  line  which  appears  as  a  noisy  curve.  The  estimated 

values  are  plotted  using  a  dashed  line  with  x's  on  it,  which  appears 

in  most  cases  as  a  solid  line  with  x's  on  it. 

3.  2  ESTIMATES  OF  BIASES  AND  SCALE  FACTORS 

The  angle  of  attack  sensor  has  a  measurement  rsuige  from 
-15*  to  45*  and  is  centered  on  a  s  i5>-  which  is  stall.  The  relationship 
between  the  true  angle  of  attack  a  and  the  measured  angle  of  attack  is 

0.  =  15*  +  (a  -  15*)  +  b  +  noise 
am  a 

where  is  the  scale  factor  and  b^  is  the  bias.  The  sideslip  angle 
sensor  has  a  measurement  range  from  -30*  to  30*  and  is  centered  on 
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/ 

.# . 


P  =  0.  The  relationship  between  the  true  sideslip  angle  /3  £ind  the 
measured  sideslip  angle  is 

j3  »  0*  +  X  a  O  -  0*)  +  bfl  +  noise 

pm  p 

The  scale  factor  X  for  the  angle  of  attack  was  estimated 
to  be  0.  78.  The  bias  b  was  estimated  to  be  2*  for  maneuvers  FlMl, 

Gw 

F2M1, . . . ,  F3M3,  and  F6M1  and  0“  for  the  maneuvers  F4M1,  F4M2, 
F4M10. 

The  scale  factor  X^  for  the  sideslip  angle  was  estimated  to 

be  0.805.  The  bias  b^  varied  from  one  maneuver  to  the  next.  For 

P 

maneuvers  FlMl,  F2M1. . . .,  F3M3  and  F6M1  it  varied  between  0. 1" 
and  1.7®.  For  the  maneuvers  in  flight  4  it  varied  from  0,8“  to  5.1®. 

Two  sensors  were  used  to  measure  the  total  airspeed. 

One  sensor  has  a  measurement  range  of  500  knots;  it's  value  is  called 
a  "coarse"  value.  The  other  sensor  has  a  range  of  100  knots  and  it's 
value  is  called  a  "fine"  value.  A  comparison  between  these  two  values 
is  presented  for  each  maneuver  in  Figures  A-1,  A-2  and  A-3  in  Appen¬ 
dix  A.  The  fine  values  are  3,  9  m/sec  larger  than  the  coarse  values 
for  maneuvers  FlMl,  F2M1, . . .,  F3M3  and  F6M1.  The  fine  and  coarse 
values  coincide  for  the  maneuvers  in  flight  4. 

The  maneuvers  exhibiting  the  3.  9  ni/ sec  difference  were  flown 
between  Nov.  5  and  Dec.  29,  1975.  The  maneuvers  of  flight  4  were 
flown  on  Aug.  26,  1976  which  is  seven  months  after  flight  3.  The  fine 
airspeed  measurement  which  has  one-fifth  the  noise  level  of  the  coarse 
airspeed  measurement  was  used  as  the  measured  variable  in  the  filter/ 
smoother.  It  has  a  bias  of  -3.9  m/sec  for  the  maneuvers  of  flights  1, 

2,  3  and  6  and  no  bias  for  the  maneuvers  of  flight  4. 
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Large  biases  were  estimated  for  the  roll  and  pitch  rates. 

The  roll  rate  bias  b  was  estimated  to  be  1065  rad/sec  which  is 

P 

-6.1  deg/sec.  The  pitch  rate  bias  b^  was  estimated  to  be  -.0479 

rad /sec  which  is  -2.  7  deg/ sec.  The  yaw  rate  bias  b^  was  estimated 

within  ±  .  02  rad/sec  which  is  within  its  large  noise  level.  An  axis 

o 

alignment  bias  of  the  rate  gyro  was  estimated  as  3  . 


The  trim  values  of  &  and  6  from  the  eighteen  maneuvers 

a  r 


were  plotted  against  each  other.  The  eighteen  points  (&  ,  6  )  lie 

&  r 


nearly  on  a  straight  line  which  crosses  the  6  -axis  at  6  »  1.8°. 

&  & 


There  is  a  cluster  of  twelve  points  at  this  crossing.  The  line  should 


cross  at  i  =0*  because  with  6  =0  trim  conditions  should  follow 
a  r 


from  6  =0  and  not  from  6  =1.8*.  This  value  of  1. 8*  was  accepted 

d  SI 


as  the  bias  on  6  .  The  other  controls  6  smd  6  had  no  biases, 
a  re 


The  bias  of  the  Fore-Aft  acceleration  a 


was  estimated 


NOSE 


to  be  0. 6  m/sec  and  the  biases  of  the  lateral  accelerations  a 


and 


c.g. 


a  are  -.  13  m/sec  .  The  bias  of  the  lateral  acceleration  a 

^NOSE  ^TAIL 

is  within  its  large  noise  level.  The  normal  acceleration  measurement  at 

2 

the  center  of  gravity  has  a  bias  of  -1.6  m/ sec  and  the  off  c.  g.  normal 


acceleration  measurements  have  a  bias  of  1.0  m/sec  . 
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SOME  OBSERVATIONS  ON  THE  T-2C  DATA 


A  correlation  analysis  of  the  data  revealed  that  the  measured 
quantities  belong  to  two  distinct  groups.  The  first  group  consisted  of 


3*  P.  c[,  9,  0,  -  and  6  .  The  second  group  consisted  of  the  re- 

9  Si 


maining  measured  quantities  (i.e.,  all  accelerations,  r,  6  ,  altitude 

r 

and  airspeed).  During  the  trim  portion  of  each  maneuver  the  noise  of 


any  member  of  the  first  group  is  80%  to  95%  correlated  with  the  noise 
of  any  other  member  of  the  first  group  and  is  0%  to  10%  correlated 
with  the  noise  of  any  member  in  the  second  group.  Similarly,  the 
members  of  the  second  group  exhibit  the  same  noise  correlation 


characteristic  with  each  other.  The  reason  that  there  are  two  groups 
and  not  one  is  because  single  channels  were  used  to  measure  the 
quantities  of  the  first  group  whereas  a  commutator  working  off  a 
single  channel  was  used  for  the  second  group.  The  reason  that  there 
is  any  group  having  such  a  high  correlation  has  not  been  explained. 

The  commutator  had  the  effect  of  adding  4  to  5  times  the  noise  level 
of  a  single  channel  level.  For  example,  p  and  r  should  normally 
exhibit  the  same  noise  level  characteristic  during  trim  since  they  have 
the  range  ±  150  deg/sec  and  since  they  are  measured  with  the  same 
type  of  instrument.  From  the  plots  of  p  and  r  in  Figure  B-1  we 
observe  that  the  noise  level  on  r  is  4  to  5  times  higher  than  it  is  on 
p.  The  variable  p  belongs  to  group  one,  and  r  belongs  to  the  second 
group.  Also,  we  note  in  Figure  B-1  that  the  noise  level  on  6^  is  con¬ 
siderably  higher  than  it  is  on  5  and  5  .  Single  channels  should  be 

6  St 

used  for  all  quantities  since  this  keeps  the  noise  levels  to  a  minimum. 
Also,  common  noise  sources  should  be  isolated  and  removed,  if  pos¬ 
sible,  from  the  instrumentation. 

The  data  exhibits  noise  from  a  nose  boom  vibrating  between 
6  and  7  hertz  with  amplitudes  as  high  as  several  degrees  in  a  and  p. 
The  off  c.  g.  acceleration  measurements  exhibit  high  vibrational  noise. 

There  are  occasional  losses  of  total  airspeed  measurement 
over  extended  intervals.  In  some  cases  the  airspeed  would  change 
abniptly  by  20  m/sec. 


4. 


IDENTIFICATION  OF  THE  T-2C  LATERAL  AERODYNAMIC 
CHARACTERISTICS 


4. 1  CONSTRUCTION  OF  SUBSPACES  FOR  LATERAL  MODELING 

The  lateral  coefficients  C  ,  C ,  and  C  are  functions  of  the 

y  *  n 

following  variables: 


OL  -  angle  of  attack 
^  -  sideslip  angle 

p  -  roll  rate 
r  -  yaw  rate 

-  aileron  deflection 

-  rudder  deflection 

s 

^  -  rate  of  sideslip 

The  states  a  and  /3  are  used  for  defining  the  subspaces  for  modeling. 

The  subspaces  in  a  are  created  by  partitioning  a  at  -4®,  -3",  -2®, . . . , 

23®,  24®,  25®,  30®,  35®  and  40®.  That  is,  one  degree  intervals  are  selec¬ 
ted  between  -4®  and  25®  and  five  degrees  intervals  are  chosen  between 
25®  and  40®.  The  eighteen  data  files  (i.  e. ,  one  for  each  of  the  maneuvers 
FlMl,  F2M1, . . . ,  F4M10,  and  F6M1)  created  by  the  estimation  phase  of 
•the  EBM  system  identification  method  are  transformed  into  thirty-two  new 
data  files  by  using  the  partitions  on  a  and  by  collecting  together  the  data 
from  all  eighteen  maneuvers  having  an  angle  of  attack  within  an  a  interval. 


The  subspaces  in  ff  are  handled  by  using  spline  functions.  Since 
the  lateral  coefficients  are  emtisymmetric  with  respect  to  the  sideslip  angle, 
a  cubic  spline  function  having  only  odd  powers  of  0  is  used.  This  ensures 
that  the  derivatives  with  respect  to  |3  are  continuous  and  are  symmetric  in 
For  example,  using  C  to  illustrate  we  have  (assuming  linear  terms 

.  y 

in  5^,  p,  r  and  jS) 
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N 


C  (a.  A  ft  =0.  p*0.  r  =  0.  ^-0)  =  C  (a).^+z:  C  (a)  •  f  (A  fl.,  .  J 
y  r  Va  y  o  i  i+'i 


^A'  i»i  V 
o  Pi 


where  and  A^  are  spline  knots  on  the  |Al  axis  satisfying 
0  s  Aj  <  Ag  and  where  f  (A,  jS^,  is  defined  as 


/O 


f  W. 


io-?/ 


A  ^  A 

'^i+l 


The  derivative  of  with  respect  to  A  is  given  by 


N 


hi 


*  iM  “y  3 ^ 

^  ^o  A- 

1 


where 


10  0  s  A  ^  Aj^ 

li  <-^>  ^i- =  li 

We  observe  from  these  derivatives  that  C  (a.  A)  Is  an  even  function 

^A 

with  respect  to  A  and  is,  therefore,  symmetric  in  A* 
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Let  3  max  (a)  denote  the  maximum  value  of  |^|  attained  at  a 
given  angle  of  attack  a  by  all  eighteen  maneuvers.  At  a  given  a  the 
subspaces  in  P  are  created  by  choosing  knots  P^,  etc.  in  the 
/3-space  [0,  P  max  (a)].  The  five  following  knot  sets  are  used  in  the 
modeling  of  the  lateral  coefficients: 


Set 

Knots  on  i0|  axis 

I 

0». 

jSmax  (a) 

n 

OS 

5S 

P  max  (a) 

in 

OS 

10» 

,  /3max  (0) 

IV 

OS 

2S 

lOS  Pmax  (a) 

V 

OS 

5S 

lOS  pToax  (0) 
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4.2  EDENTIFICATION  OP  C 

y 

As  discussed  in  Section  4. 1  the  modeling  of  is  conducted  on 
subspaces  of  the  a -space  -4®  to  40®.  The  subspaces  are  one  degree  in¬ 
tervals  [a.,  +  1®],  *  -4®,  -3®, ....  24®  and  five  degrees  intervals 

[a^,  +  5®].  =  25®,  30®,  35®. 

Let  a-^  be  the  left  hand  endpoint  of  an  a-intervsd  and  let  be 
the  right  hand  endpoint.  Let  a  be  the  midpoint: 

-1 - 

The  estimation  results  of  all  eighteen  maneuvers  having  an  angle  of 
attack  within  the  a -interval  [a^,  o^]  is  processed  using  the  stepwise 
multiple  linear  regression  (SMLR)  technique,  the  modeling  phase  of  the 
EBM  system  identification  method.  The  following  model  with  a  fixed  at 
a  ®  a  is  used  in  the  SMLR  technique  to  identify  a  model  for  the  a -interval 
<Xj]= 


(a=».  S.  P.  r,  C  W-6+C  W.®+C 


a  r 


^6  ^  ^6 


2v  y  '  2v 
■'r 


JN 

+  C  (a)  •  ^  +  C  (^  •  ^  +  E  C  (ct)  •  f  O,  P.,  P.  .)  +  other  pos- 
^P  ^  ^P  i=l  y^3  " 


sible  terms. 


The  SMLR  technique  found  the  parameters  C  (a),  C  (a),  C  (a) 

and  C  to  be  statistically  insignificant.  The  identification  of  the  remaining 
^P 

parameters  are  given  in  Table  C-1  of  Appendix  C.  For  example,  at  a  =  18.5®, 
the  identified  parameters  are: 
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C  (a  =  18. 5»)  =  .076 
^6 

r 

C  (a  *  18.5**)  *  -.67 


=  18.5“)»-4.56 

=  18.5*)  *  14.76 


where  the  knots  are  0^  -  0**,  =  5®,  and  0^  ~  18®.  Due  to  the  sparsity 

of  data  at  low  angles  of  attack  and  at  high  angles  of  attack  several 
a -intervals  were  combined  in  order  to  identify  the  parameters.  For  ex¬ 
ample,  at  the  low  end  of  the  ot-spectrum  the  following  eqtiation  form  was 
used  in  the  SMLR  program: 


•  X  (0)  +C  £)  •  X  (0) 
J  O' A  ^  l-AO  '»0\  Jff  ^  /  .«  !-«' 


(-4®,  7®)  ^0 


(-4®,  5®) 


'y 

0 

1 

'y 

0' 

1 

'y 

0 


(Sj  • 

0 

•  X 

(a)  + 

C 

(a_) 

(5®, 

6®) 

0 

(a  )  • 

’  f 

<0. 

0®, 

o 

o 

•Xu*' 

• 

X 

(a) 

1 

(-4® 

,  5®: 

(aj  ' 

-  f 

(0, 

0®, 

10®)  • 

X 

(a) 

A 

(5®. 

6®) 

{5J  • 

•  f 

(0, 

0®. 

12®)  • 

X 

(a) 

0 

(6®. 

7®) 

(6®,  7®) 


where  x  (<^)  is  the  characteristic  function  of  the  set  [a.  b]  defined  as 


(a, 

X  (0.)  =  ^ 

(a,  b) 


1  if  a  i  as  b 


0  if  a<a  or  b<a 
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and  where  =  0.5%  -  5.5“  and  =  6.5“.  The  estimate  of  a 

parameter  holds  over  the  range  of  the  characteristic  function  that 
multiplies  it.  The  value  of  C  was  identified  as  .118.  Since  the 

range  of  the  characteristic  funcfion  that  multiplies  C  is  from 

-4“  to  7“,  the  estimate  of  C  -  .  118  is  tabulated  for  each  of  the 

intervals  f-4,  5“L  [5“,  6®]  anJ"  [6®,  7®]  listed  in  the  first  column  of 
Table  C-1.  The  other  six  parameters  are  tabulated  in  Table  C->1 
under  the  column  headings  of  C  and  C  .  Characteristic  ftmctions 


/s: 


were  used  also  at  the  high  angles  of  attack  to  obtain  the  identified  param¬ 
eters  of  Table  C-1  over  the  a -range  between  24®  and  40®. 

The  aerodynamic  derivative  C  (a)  which  is  given  in  Table  C-1 

is  plotted  in  Figure  4. 2-1  where  it  is  compared  with  the  wind  tunnel  model. 


model  and  the  wind  tunnel  moSel 


The  parameters  given  in  Table  C-1  were  used  to  plot  Figures 
4. 2-2  through  4.2-6.  These  figures  compare  the  identified  models  of 
C  (a)  and  of  C  (a)  with  the  wind  tunnel  model  for  6=0  and  for 

y  y  a  ^ 

13  =  0*,  1®,  ...,10"  The  derivative  C^,  08)  for  a  =10.5®,  11.5®,..., 
23.5®,  24.5®,  27.5®  is  presented  in  Fi^ires  4.  2-5  and  4.  2-6.  Observe 
in  Figure  4. 2-6  the  match  between  the  identified  C  O  and  the  wind 
tunnel  for  a  =  20.5®,  21.5®  and  22.5®  as  well  as  for  a  =  27.5®. 


The  SMLR  technique  estimates  in  addition  to  the  parameters  of 

Table  C-1  a  value  for  Cy^  (5)  at  each  of  the  a  -  intervals.  The 

parameter  C  (a)  is  the  value  of  the  side  force  coefficient  C  when¬ 
ce  —  ^ 

ever  all  states  and  all  controls  are  zero  except  a  =  a.  The  estimated 

—  -4 

values  of  C  (a)  were  of  the  order  10  which  being  negligible  were 
taken  as  zero. 


'  {  _  I  ni.PH«|tocqi  I  I  I 

Fig.  4.2-2  Comparison  between  the  C  and  the  C 

^  y  ( 

wind  tunnel  model  for  0  >  0**.  1°,  2”,  3°. 


;  flUPHfitOEO)  ; 
identified  models'  and  the 


Fig.  4.2-4  Comparison  between  the  C^,  and  the  identified  models  and  the  wind 

tunnel  model  for  /3  =  8*.  9*,  10*  ^ 
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4.3  IDENTIFICATION  OF 

The  modeling  of  is  conducted  on  subspaces  of  the  a -space 
-4®  to  40®.  The  subspaces  are  one  degree  intervals  =  1®], 

=  -4®,  -3®, . . . ,  24®  and  five  degrees  intervals  [ot,  a.  +  5®],  =  25®, 

30®,  35®.  Let  a  denote  the  midpoint  of  an  a-intervsd..  The  following 
model  with  a  fixed  at  a  =  a  is  used  in  the  SMLR  program  to  identify 
the  aerodynamic  coefficient  for  the  ot  -interval  being  modeled: 

C^(5:  .  jS,  p,  r,  6^,  h  =  (5)  .  6^  +  (5)  •  6^ 

6  6 
a  r 

p  r  P 

_  N  _ 

+  C.  (a)  •  ^  L  C.  (a)  •  f  03,  /3.,  +  other  possible  terms. 

i=l  \2  ^  ^  ^ 

o  ^ 

The  parameter  C, ,  (a)  was  determined  by  the  SMLR  technique  to  be 

P 

statistically  insignificant.  The  identified  values  of  the  other  parameters 
are  presented  in  Tables  D-1  and  D-2  of  Appendix  D.  For  example,  at 
a  =  14.5®,  the  identified  parameters  are: 

C.  (a  =  14.5®)  *  -.148  per  rad 
a  _ 

C.  (01  *  14.5®)  =  .0127  per  rad 
r__ 

(a*  14.5®)  =  -.225  per  rad 
P 

(5  *  14.5®)  =  .310  per  rad 
r 

C,  (a  »  14.5®)  =  -.037  per  rad 

\ 
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_  (a  =  14.5®)  =  14.23  per  rad^ 

C,  (a  =  14.5®)  =  -11.21  per  rad^ 

2  -  3 

C.  (a  =  14.5®)  =  102.2  per  rad 

^  3 
jS 

where  the  knots  are  -  0®,  ^  =  2®,  ^  =  10®  and  ]3  *  15®. 

X  A  u  4 

Due  to  the  sparsity  of  data  at  low  and  high  angles  of  attack, 
characteristic  functions  (as  discussed  in  Section  4.  2)  were  used  to 
combine  several  a -intervals  in  the  modeling  of  For  instance, 
the  a -intervals  below  a  =  6®  were  combined  using  characteristic 
functions  as  weU  as  those  above  a  =  22®.  In  addition  to  using  charact¬ 
eristic  functions,  the  following  procedure  was  used  to  identify  the 

slopes  of  (O')  and  C.  (o)  at  low  ^d  high  angles  of  attack:  we  use 
1^.  r 

C  (a)  .  to  itluBtrate  it.  The  value  of  C  (a  »  21.5®)  has  been  identi- 

p  p 

fled  by  processing  the  data  in  the  a -interval  [21®,  22®].  It  is  desired 

to  identify  the  slope  of  C.  (i)  for  a  in  the  range  22®  to  40®.  This 

is  accomplished  by  modifying  the  term  C .  (a)  •  2^  as 

P 

[C  (a  =  21.5®)  +  A-^2L^^^]*  ^  •  X(a) 

^p  2V  ^22®,  40») 

and  identifying  the  slope  A  by  processing  the  data  in  all  a>intervals 

between  22®  and  40®.  The  slope  A  is  then  used  to  evaluate  (a)  for 

P 

any  a  in  the  range  22®  to  40®; 


Cjj  (a)  =  (a*  21.5®)  + A  -  22®  s  as  40® 
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• 


This  procedure  was  used  in  obtaining  the  tabulated  values  in  Table  D-2 
for  a*22.5».  23. 5».  24. 27. 32. 5“  and  37.  5»  and  in  obtaining 

(a)  for  a  in  the  range  -4®  to  10®. 
r 

The  identified  dynamic  derivative  (a)  which  is  tabulated  in 
Table  D-2  is  plotted  in  Figtire  4.3-1  where  i^^is  compared  with  the 
theoretical  prediction  model  of  Bihrle  Applied  Research,  Ref.  [2].  At 
low  a  the  identified  model  of  (a)  has  sUghtly  smaller  values  than 


.87  18.00  .  Z 

vftUPHfl:0EQ) 


Fig.  4.3-1  Comparison  between  the 


and  the  theoretical  prediction  model  of  Bihrle 
Applied  Research 


the  theoretical  prediction  model.  The  two  curves  cross  at  about  15.5® 
and  the  identified  model  hovers  about  the  zero  axis  between  16®  and  20® 
whereas  the  theoretical  prediction  goes  positive  and  reaches  a  peak  of 
0.30  per  rad.  The  two  curves  agree  for  a  above  20®. 
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and  the  theoretical  prediction  model  of  Bihrle 
Applied  Research 


steeper  slope  at  low  angles  of  attack  but  drops  off  more  sharply  at 
stall  than  the  theoretical  prediction  model.  Above  a  =  18”  the  identi¬ 
fied  model  lies  slightly  imder  the  theoretical  prediction. 

The  tabulated  values  of  C,  (a)  in  Table  D-1  are  plotted  in 

a 

Figure  4.3-3,  together  with  the  wind  tunnel  model  for  j3  =  0®.  The 
identified  model  of  C,  (<x)  compares  well  with  the  wind  hmnel  model 

a 

at  low  a.  Its  drop  in  controllability  above  stall  coincides  with  the  wind 
tunnel  model's  drop.  Between  18®  and  26®  angles  of  attack  the  identified 


r.zia-: 


I 


Fig.  4.3-3  Comparison  between  the  C^  identified 

a 

model  and  the  wind  tunnel  model 


model  has  a  slightly  smaller  magnitude  than  the  wind  tunnel,  but  they 
coincide  above  30”. 

The  identified  values  of  C.  (a)  in  Table  D-1  are  plotted  in 


Figure  4.3-4  where  it  is  compared  with  the  wind  tunnel  model  for 
P  =  0®.  They  compare  well  except  for  a  between  stall  and  20”  where 
the  wind  tunnel  model  seesaws  to  double  its  magnitude  and  the  identi¬ 
fied  model  seesaws  about  zero. 


Fig.  4.3-4  Comparison  between  me  C.  identifiec 

6 

r 

model  and  the  wind  tunnel  model 


The  other  identified  parameters  of  Table  D-1  were  used  to 
plot  Figures  4.  3-5  through  4. 3-9  where  C  (a)  and  C.  (a),  -4®  s  as  40® 

P 

are  compared  with  the  wind  tunnel  model  for  /3  =  0®,  1®, . . . ,  lO®  and 

where  (P).  0  s  ^  s  10®,  is  compared  with  the  wind  tunnel  model  for 
P 

a  10.5®,  11.5®,...,  23.5®,  24.5®,  27.5®. 

The  SMLR  technique  estimates  in  addition  to  the  parameters  of 

Tables  D-1  and  D-2  a  value  for  W  at  end  of  the  a-intervals.  The 

o 

parameter  (a)  is  the  value  of  the  rolling  moment  coefficient  C 
o 

whenever  all  states  and  all  controls  are  zero  except  o^»a  =  ^  The 
estimated  values  of  were  of  the  order  10 to  10 which  being 

negligible  were  taken  as  zero. 


.33 '  at.oo 
8ETR(0E0I 


I  8ETBIOEOJ  ;  '  .  _1.  KTMOEOI  i 

4.  3-8  Comparison  between  the  identified  model  and  the  wind  tunnel  model  for 
a«10.5*,  11.5",  12. 5».  13. 14. 5«,  15. 5»,  16. 5»,  17.5“ 
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4.4  IDENTIFICATION  OF  C 

n 

The  modeling  of  is  conducted  on  subspaces  of  the  o^space 
-4“  to  40®.  The  subspaces  are  one  degree  intervals  [ot.,  +  1*], 

=  -4®,  -3®, . . . ,  24®  and  five  degrees  intervals  [a^, 

=  25®,  30®,  35®.  Let  a  denote  the  midpoint  of  an  a-interval.  The 
following  model  with  a  fhced  at  a  =  a  is  used  in  the  SMLR  program  to 
identify  the  aerodynamic  coefficient  over  a  given  o-inteirval: 


(a,  fi.  p,  r,  6^,  fi>  =  C  (a)  .  6^  +  (a>  •  6^ 

0  0 
a  r 


+  c  (i)  •  +  c_  (5)  •  ^  +  C_.  (a)  . 


2V 


2V 


^18 


2V 


N 

+  c  (S)  •  j3+  E 
\  i=l 


f  (jSt  +  other  possible  terms. 


The  parameter  C  (a)  and  C  (a)  were  determined  by  the  SMLR 

program  to  be  statistically  insignificant.  The  identified  values  of  the 
other  parameters  are  presented  in  Tables  E-1  and  E-2  of  Appendix  E. 
For  example,  at  a  =  13.5®,  the  identified  parameters  are: 


C  (a  =  13.5®) 

r 

C  (5  =  13.5®) 

“p 

C  (a  =13,5) 
n 

r 

C  (a)  =13.5®)  = 

C  (a  =  13.5®)  = 


-.032  per  rad 
-.098  per  rad 
-.  1816  per  rad 
.  052  per  rad 
4.18  per  rad^ 
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C  (d  s  13.5®)  =  -15.93  per  rad^ 

^  3 

4 

C  (£  =>  13.  S')  •  32. 19  per  rad* 

n  _ 

B 

where  the  knots  are  0^^  -  0®,  0^  ~  5®,  0^  -  10®  and  0^  =  16®. 


Due  to  the  sparsity  of  data  at  low  and  high  angles  of  attack. 

characteristic  functions  (as  discussed  in  Section  4.2)  were  used  to 

combine  several  a>intervals  in  the  modeling  of  C^.  In  addition  to 

using  characteristic  functions,  the  procedure  for  estimating  slopes 

described  in  Section  4.3  was  used  to  identify  C  (oi),  C  (a)  and 

n  n 

P  r 

C  (a)  for  a  above  25®.  The  procedure  was  used  to  identify  (a) 

and  C  (a)  over  the  a  range  25®  to  40®  and  C  (a)  over  the  a  range 
n,  n 

6  r 

r 

30®  to  40®. 

The  identified  dynamic  derivative  (a)  which  is  tabulated  in 
Table  E-2  is  plotted  in  Figure  4.4-1  where  i^’is  compared  with  the 
theoretical  prediction  model  of  Bihrle  Applied  Research,  Ref.  [2].  The 
identified  model  agrees  with  the  theoretical  prediction  at  angles  of  attack 
below  stall.  The  identified  model  reaches  its  highest  positive  value  of 
.113  at  a  =  16.5®  whereas  the  theoretical  prediction  model  continues 
to  grow  reaching  its  peak  of  0. 365  at  a  =  18®.  After  reaching  a  peak  at 
a  s  16.5®  the  identified  model  drops  off  to  a  value  of  -.076  at  a  =  20.5®, 
rises  to  .048  at  a  *  22.5;  then  drops  to  -.051  at  a  =  23.5  before  level¬ 
ing  off  about  -.075  over  the  at  range  [25®,  40®].  The  theoretical  predic¬ 
tion  model,  on  the  other  hand,  drops  from  its  peak  at  a  =  18®  to  .  12  at 
a  a  20.5®,  climbs  to  .  14  at  a  =  22®,  drops  slightly  to  .  13  at  a  »  24®,  then 
rises  to  .  195  at  a  =  30®  before  dropping  to  .  06  at  a  =  38®.  A  significant 
difference  between  these  two  models  is  that  over  the  a  range  [20®,  35®] 
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Fig.  4.4-1  Comparison  between  the  identified 

model  and  the  theoretical  pridiction  model 

of  Bihrle  Applied  Research 

the  theoretical  prediction  model  has  values  greater  than  .  10  while  the 
identified  model  has  values  around  -.075  except  for  the  sharp  rise  and 
decline  between  a  =  21”  and  23”  where  it  peaks  at  .048.  This  difference 
is  highly  significant  as  has  been  demonstrated  by  the  findings  of  the 
simulation  study  of  [1].  Therein,  the  initial  conditions  and  the  recorded 
control  time  histories  of  the  actual  T-2C  flight  data  were  used  to  drive 
the  synthetic  T-2C  model  composed  of  the  wind  tunnel  model  and  the 
theoretical  prediction  model  of  Bihrle  Applied  Research.  The  synthetic 
responses  were  drastically  different  from  the  actual  flight  data  until  a 
change  in  the  theoretical  prediction  model  of  (i.  e. ,  the  same  as 
given  in  Figxire  4.4-1)  was  made.  For  example.^one  synthetic  maneuver 
reached  60”  in  o-  while  the  real  data  did  not  go  beyond  a  =  25”.  Other 
synthetic  maneuvers  exhibited  large  yaw  and  roll  spin  rates.  An  inspec¬ 
tion  of  the  states  revealed  a  greater  degree  of  instability  for  the  theoretical 


prediction  model  of  .  By  trial  and  error  a  constant  value  of 

C  *  06  was  found  t<P give  stability  to  the  simulated  data  and  to 

n 

provide  data  which  agreed  better  with  the  actual  flight  data.  Section 
2  of  Ref.  [1]  contains  a  more  detailed  discussion  of  the  problem. 
Therein,  note  also  the  results  in  Figure  4.17.  Therefore,  the  identi¬ 
fied  model  of  obtained  by  processing  the  actual  T-2C  flight  data 

agrees  with  the  nSodification  of  obtained  by  trial  and  error  in 

Ref.  [1].  ^ 

The  tabulated  values  of  C  (a)  in  Table  E-2  are  plotted  in 

n 

Figure  4.4-2  where  it  is  comparetf with  the  theoretical  prediction 
model  of  Ref.  [2].  Amazingly,  the  theoretical  prediction  model  agrees 
well  with  identified  model.  There  appears,  however,  to  be  a  difference 
between  the  curves  of  about  1®  shift  in  a. 


/ 


The  tabulated  values  of  C  (  a)  in  Table  E-1  are  plotted  in 

Figure  4.4-3  where  it  is  compared  with  the  wind  tunnel  model  for 
/3  =  0°.  The  identified  model  indicates  a  slightly  smaller  performance 
than  that  expected  by  the  wind  tunnel  results. 
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Fig.  4.4-3  Comparison  between  the  identified 


model  and  the  wind  tunnel  model 


The  other  identified  parameters  of  Table  E-1  were  used  to  plot 

Figures  4.4-4  through  4.4-8  where  C  (a)  and  C  (a),  -4®  s  a  s  40®, 

n  ng 

are  compared  with  the  wind  tunnel  model  for  -  0®,  1®,  — ,  10®  and 
where  C  (|3),  0  s  ^  s  10®,  is  compared  with  the  wind  tunnel  model  for 

a  =  10.5®,  11.5®,.,.,  23.5®,  24.5®,  27.5®.  The  rudder  control  equals 
zero  in  these  figures. 


The  SMLR  technique  estimates  in  addition  to  the  parameters 

of  Tables  E-1  and  E -2  a  value  for  C  (5)  at  each  of  the  a -intervals. 

n 

o 
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The  parameter  C  (a)  is  the  value  of  the  yawing  moment  coefficient  C 
o  “ 

whenever  all  states  and  controls  are  zero  with  the  exception  of  a  -a. 

—  -4  _5 

The  estimated  values  of  C  (a)  were  of  the  order  10  to  10  which 

n 

being  negligible  were  taken  as  zero.  A  linear  term  involving  p  was 
identified  to  have  a  nearly  constant  value  of  1. 1  for  all  a- intervals. 
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5.  IDENTIFICATION  OF  THE  T-2C  LONGITUDINAL 
AERODYNAMIC  CHARACTERISTICS 


5. 1  CONSTRUCTION  OF  SUBSPACES  FOR  LONGITUDINAL 
MODELING 


The  longitudinal  coefficients  C  ,  C  , 

X  z 

of  the  following  variables: 


and  C  are  functions 
m 


a  •  angle  of  attack 

0  -  sideslip  angle 

q  -  pitch  rate 

6  -  elevator  deflection 

e 

a  -  CL  rate 

The  states  a.  0  and  are  used  for  defining  the  subspaces  for  model¬ 
ing.  Subspace  modeling  is  conducted  in  two  steps  for  the  modeling  of  C  . 

m 

The  first  uses  a  partitioning  of  the  a  space  [-4®,  40®]  together  with  spline 
functions  to  handle  the  subspaces  of  the  0  and  the  6  space.  The  sub- 
spaces  in  a  are  generated  by  partitioning  the  a  space  into  one  degree 
intervals  +  1®],  *  '4®,  -3®, . . . ,  24®  suid  five  degrees  intervals 

[Ok.  +  5®],  =  25®,  30®,  35®.  The  eighteen  estimation  data  files  (i.e., 

one  for  each  of  the  maneuvers  FlMl,  F2M1, . . .,  F4M10  and.FBMl)  gen¬ 
erated  by  die  estimation  phase  of  the  EBM  system  identification  method 
are  transformed  into  thirty-two  new  estimation  data  files  for  modeling  by 
using  the  partitions  on  a  and  by  collecting  together  the  data  from  all 
eighteen  maneuvers  having  an  angle  of  attack  within  a  given  a-interval. 

In  this  first  step  the  subspaces  in  0  and  6^  are  handled  by  using  spline 
functions.  Since  a  longitudinal  coefficient  is  symmetric  with  respect 
to  the  sideslip  angle,  a  quadratic  spline  fimction  as  a  function  of  |/9|  is 
used. 
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The  second  step  uses  a  partitioning  of  the  6  space  [-28",  2"] 

together  with  spline  functions  to  handle  the  subspaces  of  a .  The  space 

of  6  is  partitioned  into  one  degree  intervals  [6,6  +  1"],  6  =  -28' 

e  e^  e^^  e^ 

-27", . . . ,  1".  The  eighteen  estimation  data  files  generated  by  the  esti¬ 
mation  phase  of  the  EBM  method  are  transformed  into  thirty  new  esti¬ 
mation  data  files  for  modeling  by  using  the  partitions  on  6^  and  by 
collecting  together  the  data  from  all  eighteen  maneuvers  having  an 
elevator  control  deflection  within  a  given  6^ -interval.  In  this  second 
step  the  subspaces  in  a  are  handled  by  using  spline  functions. 
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5.2  IDENTIFICATION  OF  C 

m 


The  estimated  values  of  C  (t)  are  compared  with  the  wind 

m 

tunnel  values  in  Figures  F-1  through  F-14  of  Appendix  F.  The  elevator 

control  6^  was  used  to  subdivide  die  C^(t)  estimates  of  all  sixteen 

maneuvers  into  data  subsets.  The  first  data  subset  consists  of  the  C  (t) 

m 

estimate  for  which  (t)  is  a  value  between  -1.0**  and  0.0*.  The  second 

e 

data  subset  is  over  the  elevator  control  interval  of  -2.0“  to  -1. 0“.  The 
data  subsets  cover  the  entire  range  of  the  elevator  control  in  one  degree 
intervals  with  the  last  data  subset  covering  the  interval  from  -28.0“  to 
-27.0“.  Each  data  subset  was  ordered  with  respect  to  the  angle  of  attack 
variable  and  then  plotted  as  a  figure  in  Appendix  F.  The  wind  :tnnnel 
values  were  generated  by  evaluating  the  wind  tunnel  model  using  the  esti¬ 
mated  states  and  controls.  Note  that  the  range  of  the  angle  of  attack 
depends  on  the  elevator  control  interval.  For  example,  for  the  interval 

of  6  »  -1*  to  0“,  the  range  on  a  is  from  -1.36“  to  18. 94“  whereas  for 

© 

the  interval  of  6^  =  -27“  to  -26“  the  range  of  a  is  from  15.51“  to  37.54“. 

This  is  to  be  expected  since  during  maneuvers  a  small  elevator  control 

setting  tends  to  produce  low  alpha  and  a  large  negative  setting  generates 

high  alpha.  The  curves  have  dependency  on  all  states  by  nature  of  the 

flight  maneuvers.  The  fair  agreement  between  the  real  data  values  and 

the  wind  tunnel  values  indicates  a  likeness  between  the  actual  aircraft 

and  its  wind  tunnel  model.  There  is,  however,  a  noticeable  diHerence. 

At  angles  of  attack  below  stall,  the  real  data  values  azte  above  the  wind 

tunnel  values  over  the  elevator  control  range  -4“  to  0“;  they  coincide 

over  »  -6“  to  -4“  auid  they  lie  below  the  wind  tunnel  values  over 

6  =  -28®  to  -7“.  The  models  that  we  identify  for  should  exhibit 
e  m 

this  same  phenomenon. 


Oir  objective  in  modeling  is  to  identify  a  model  of 

for  each  of  the  6  -subspaces  [6,6+  1“],  6  =  -28“,  -27“ 

e  e.  e.  e. 


$  •  •  •  $ 


That  is,  we  model  "the  real  data  values"  curves  given  in  Appendix  F. 

Let  6  be  the  midpoint  of  a  6  -interval.  An  equation  for  modeling 
©  ^ 

C  on  this  6  -interval,  valid  at  6  =6  is: 
me  e  e 

C  «!.  /S.  q.  4  6  "  C  (J  )  c  m  ■  il. 

o  q  2  V 

3 

+  C„  |«|  +  2  (<i»  •  h  0,  ?  ,  ?  ,) 

“s  j-1  “,2  ’ 

O  p. 


^  «e>  •  ^  ^  <=0.  '‘e>  •  8  '‘iH-l* 

%  - 


where  0.,  0.,  /5,,  are  the  knots  on  the  l/3l -axis  satisfsdng 
1  2  o  4 

0  S  A  0o"^^A»  where  h  (0,  0.,  0.  .)  is  a  quadratic  spline 

X  2  u  4  3  3^X 

function  defined  as 


h  tf,  S..  ■ 


i  (1^1  -  fy 


19| 


where  a,,  <x  . . . ,  a  are  the  knots  on  the  a-axis  satisfying 

-4"  s  a  j  <  • . .  <0-^,  where  g  (a,  a.,  is  a  quadratic  spline 

function  defined  as 


a<a^ 

1 


g  «i,  a^, 


2  ,  >2 


a  fi  asa. . , 
i  x+1 
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The  parameters  C  (a),  C  (a),  C  (a),  C  (a)  in  the 
m  m*  m-  m  _ 

^  Pq 

j 

preceding  equation  are  functions  of  a  and  not  6  .  Thus,  these  param- 
eters  can  be  better  estimated  by  processing  the  estimation  data  partitioned 
in  a -intervals.  Let  a  be  the  midpoint  of  an  cl -interval.  An  equation 
for  modeling  on  this  a -interval,  valid  at  a  *"5.  is: 

(a  s  a,  q,  a,  )  =  C  &)+C  (5)  •  q  +  C  £)  •  olc 

m  e  m  m  m»  —  — 


+  C  &  •  1/31  +  E  C  (a)  •  h  (^,  /3  ,  /3  ) 

m-  m  3  3+1 


+  ^  •  F  (fi  .  .  6  ) 

6  1  k=l  .2  ®  k  k+1 

o  ej^ 

where  6  ,6,6  ,6  are  the  knots  on  the  6  -axis  satisfying 

®1  ®2  ^ 

-28®  S6  <6  <6  <6  and  where  F  (6  ,  6  ,6  )  is  a  quadratic 

®1  ®2  ®3  ®4  ®  ®k  ®k+l 


spline  function  defined  as 


V’Vl’  ii“e-V' 


6  <  6 
e  e. 


6  5  6  s  6 

®k  ®  ®k+l 


,(6  -6)  (6  -46  -^5  )  6  <6 

e,  . ,  e,  e  2  e,  2  e,  , ,  e.  e 

k+1  k  k  k+1  k+1 


The  SMLR  technique  has  been  used  to  process  the  data  in  each  a -interval. 

The  identification  of  !he  parameters 

(a),  (a),  C  W,  3  =  1,  2,  3,  are  given  in  Tables  1-2 

m  m.  m-  m  „ 

q  ^  ^o  ^ 

and  1-3.  Characteristic  functions  which  are  discussed  in  Section  4.  2  used 


used  in  identif3dng  these  parameters.  At  a  =  14.  5°  we  note  in  these 
tables  that 

C  (a  ®  14.5®)  *  3.97  per  rad 
m 

q 

C  (0.  a  14.5®)  =  -.76  per  rad 

m. 

C  (a  *  14.5®)  ~  .66  per  rad 


C  (a  a  14. 5®)  a  -12.  91  per  rad' 
“  2 

(a)  a  14.5®)  a  7.98  per  rad^ 


i 


C  (a  a  14.5®)  — 17.34  per  rad 
m  ^ 


where  -  0®,  ^3  ®  ^4  “  • 

The  parameters  (a)  and  C  (a),  k  a  i  2,  3  were  used  to  plot 

“^6  “6 
®o  ®k 

in  Figures  5.2-14  suid  5.2-15  the  derivative  C  (6  )  for  the  midpoints 

m  e  ^ 

0 

_  e 

a  a  9.5®,  10.5®,  ....  24.5® 

Using  the  SMLR  technique  and  the  estimates  of  the  parameters 

given  in  Tables  1-2  and  1-3,  the  estimation  data  in  each  6 ^-interval  are 

processed  to  estimate  the  derivatives  C  (6  ),  C  (6  ),  C  (6  ), 

mo  e  m^  e  m.  ^  e 

o 

i  a  1,  2, . . . ,  6.  The  results  are  given  in  Table  I-l.  For  example,  at 

6  a  -17,5®, 
e 

C  (i  =  -17.5®)  a  .204 
m  e 
o 


C  (6  =  -17.5®)  a  -1.18  per  rad 
m  e 
a 

o 
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C  (6  =  -17.5“)  *  9.09  per  rad 

IZl  ® 

C  (fi  *  -17.5“)  *  77.8  per  rad^ 
m  ,  e 

"2. 

C  (6  =  -17.5*)  *  -11.54  per  rad 
m  -  e 

cC 

3  -  2 

C  (6  =  -17. 5*)  *  36.  96  per  rad 
m  2  e 
a 

4 

C  (6  *  -17.5*)  =  -2.  63  per  rad^ 
m  „  e 

“5_ 

C  <6  *  -17. 5*)  *  17.  72  per  rad^ 
m2  e 

CL 

6 

The  identified  model  of  C  is,  therefore, 

nx 

(a,  $.  q,  a,  &  )  *  (6  )  +  C  (a)  •  q  +  C  (a)  •  a  c 

m  e  m  e  m  m  •  - 

o  q  2  V 

*  c  IfL)  .  1^1  E  C  (a)  .  h  0.  ?  ?  ) 

Po  (S?  ‘ 

6 

+  C  (6  )  •  (a-a  )+  E  C  (6  )  •  g  (0,  Ck,  a  ) 
m^  e  1  m2  e  i  i+l 

O  CL. 

I 

where  the  parameters  C  (6),  C  (6),  C  (6),  i  =  1,  2,...,  6 

me  m^  e  m  _  e 
o  a  „2 

o  ct 

are  computed  by  making  interpolations  in  Table  I-l  and  where  the 

parameters  C  (a),  C  (a),  C  (a)  and  C  (a)  are  calculated  by 
m  m»  m-  m  _ 

1  ^ 

] 

making  interpolations  in  Tables  1-2  and  1-3. 


The  dynamic  derivative  (a)  which  is  given  in  Table  1-3  is 
plotted  in  Figure  5.2-1  where  it  is  coinpared  with  the  theoretical  predic¬ 
tion  model  of  Ref.  [2].  The  identified  model  of  C  varies  between 

m 

q 
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-11.4  and  -8.0  at  a  below  4.5®.  It  then  climbs  slowly  to  its  peak  of 
5.  61  at  a  =  13. 5®  before  falling  to  -7.  7  at  a  =  18. 5  where  it  levels  off 
to  values  between  -7.  7  and  -6. 3  for  a  beyond  18.  5®.  The  theoretical 
model  holds  almost  level  with  a  value  of  -10. 3  below  a  =  14®.  It  climbs 
sharply  to  a  peak  of  6 . 9  at  a  =  1 6®  before  dropping  slowly  to  -4.5  at 
a  =  30®. 


q 

model  and  the  theoretical  prediction  of 
Bihrle  Applied  Research. 


The  identified  model  has  positive  values  between  11®  and  16®  alpha 
whereas  the  theoretical  model  has  positive  values  between  15.5®  and 
22. 5®  alpha. 

The  dynamic  derivative  C  (a)  given  in  Table  1-3  is  plotted 

m » 

in  Figure  5.2-2  where  it  is  compare  with  the  theoretical  prediction 

model  of  Ref.  [1].  The  identified  model  of  C  is  a  constant  of  -.23 

ni» 

below  a  =  13.5®  where  it  drops  off  to  -2.  58  at  a  =  17. 5®.  It  then  rises 


with  the  same  slope  as  the  theoretical  model  until  a  «  20.5®  before 
tapering  off  to  -1. 8  at  a  =  40®.  The  theoretical  model  of  C  is 
nearly  a  constant  of  -3.3  below  ct  =  16®,  and  it  is  zero  above  =  24®. 


:  >ncMw  gartf*  MUttMiat 


!  •  -r  ■  -  i 


.oa  I  31.33 


10.87  IB.OttJ 

i  flliPHfllOEQ) 


Fig.  5.2-2  Comparison  between  the' 


identified 


model  and  the  theoretical  prediction  of 
Bihrle  Applied  Research 


The  identified  parameters  of  Table  I-l  were  used  to  plot 

Figures  5.2-3  through  5.2-9  where  C^(a)  and  C  (a)  are  compared 

™a 

with  the  wind  tunnel  model  for  6^  between  -27.5®  and  -0.5®.  The  identi¬ 
fied  parameters  of  Table  1-2  were  used  to  plot  Figures  5.  2-10  through 
5. 2-12  where  (0)  are  compared  with  the  wind  tunnel  model  for  a 
between  1.5®  and  24.5®  . 

On  the  left  hand  side  in  Figure  5.  2-3,  6^  varies  top  to  bottom 
from  6^  a  -7. 5"  to  -4. 5®;  on  the  right  hand  side  it  varies  top  to  bottom 
from  6^  •  -3.5®  to  -0.5®.  Note  as  6^  decreases  from  0.5®  to  -6.5® 


that  the  gap  between  the  identified  and  the  wind  tunnel  model  re¬ 
duces  to  nearly  zero.  Also  note  that  the  identified  model  is  above  the 
wind  timnel  model.  Then  in  Figure  5.  2-4  we  note  that  as  6^  decreases 
further,  the  wind  tunnel  model  is  above  the  identified  model  and  it  opens 
up  a  gap  for  a  below  stall.  In  Figure  5. 2-6  we  see  that  the  identified  . 
model  of  C  (a)  compares  well  with  the  wind  tunnel  model  for  a  above 

stall  and  for  6  between  -27.5“  and  -24.5“. 
e 
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.17 


RLPHRIOEO) 


Fig.  5.2-5  Comparison  between  the 


6g  *  -23.5*,  -22.5",  -21.5“,  -20.5% 


^ _  j  RLPHnfOEO)  i 

identified  model  and  the  wind  tunnel  model  for 
20.5%  -19.5“,  -18.5“,  -17.5“,  -16.5“ 


Fig.  5.2-8  Compariaon  between  the  identified  model  and  the  wind  tunnel  model  1 

6  --15.5*.  -14.5*.  -13.5»,'=^-12.5».  -11.5%  -10.5%  -9.5%  -8.5% 

e 
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Fig.  5. '2 -9  CompariBon  between  the  identified  model  and  the  wind  tunnel  model  for 

*  -23.5%  -22.5%  -21.5%  °^-20.5%  -19.5%  -18.5%  -17.5%  -16.5" 
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5.3 


RESULTS  ON  THE  C  MODEL 

z 

The  estimated  values  of  C  (t)  are  conmared  with  the  wind 

z 

tunnel  values  in  Figures  G-1  through  G*'14  in  Appendix  G.  The  ele¬ 
vator  control  5  was  used  to  subdivide  the  C  (t)  estimates  obtained 
e  z 

from  the  eighteen  maneuvers  into  data  subsets.  These  subsets  were 
ordered  with  respect  to  alpha  and  then  {dotted  in  the  figures  of  Appen¬ 
dix  F,  labeled  as  real  data  values.  The  figures  show  a  comparison 
with  the  wind  tunnel  nsodel.  The  data  in  these  subsets  were  never 
processed  to  generate  an  equation  model  of  C  >  but  the  figures  show 
that  the  model  of  the  actual  aircrsdt  is  close  to  that  given  by  the  wind 
tunnel  analysis.  In  these  figures  observe  that,  below  stall,  the  real 
data  values  almost  coincide  with  the  wind  tunnel  values  which  implies 

that  the  actual  C  (a)  is  close  to  the  wind  tunnel  model.  Aroimd  stall 
CL 

there  appears  to  be  a  hysteresis  effect  in  the  real  data  that  is  not  ac¬ 
counted  for  by  the  wind  tunnel  model.  For  example,  this  can  be  observed 
in  Figvire  G-8  for  6^  =  -14.5*, 

The  curves  given  in  Appendix  G  constitute  the  extent  to  which 

we  compared  the  wind  tunnel  model  with  that  which  arises  out  of  the 

real  data  for  the  coefficient  C  .  Time  permitting,  we  could  have 

z 

identified  an  equation  model  for  C  by  using  SMLR  technique  just  as 

z 

we  did  for  C  . 

m 
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5.4  RESULTS  ON  THE  C  MODEL 

X 

The  estimated  values  of  C^(t)  are  compared  with  the  wind 
tunnel  values  in  Figures  H-1  through  H-14.  In  these  figures  the 
aerodynamic  coefficient  is  shown  as  a  function  of  alpha  for 
Hxed  one  degree  intervals  of  the  elevator  control.  The  range  of 
the  elevator  control  covers  -28*  to  0*.  The  figures  were  obtained 
as  follows:  A  one-degree  £^-interval  was  selected.  All  eighteen 
maneuvers  were  searched  to  Hnd  all  C^(t)  estimates  having  a  cor¬ 
responding  6^  within  the  selected  6^ -interval.  The  C^(t)  estimates 
were  then  ordered  according  to  alpha  and  plotted  as  a  "real  data  values" 
curve  in  Appendix  H.  The  estimated  states  were  used  to  evaluate  the 
wind  tunnel  model  in  order  to  get  a  comparison.  These  curves  have  a 
dependency  on  q  and  a  built  into  them  by  the  nature  of  the  maneu¬ 
vers. 

The  curves  in  Appendix  H  exhibit  some  differences  between 
the  wind  tunnel  model  and  the  actual  aircraft.  An  hysteresis  effect 
can  be  observed  in  the  figures,  especially  Figure  H-8.  The  real  data 
values  are  above  the  wind  tunnel  values  at  low  alpha  and  at  high  alpha. 
Near  stall  they  are  below  the  wind  tunnel  values. 

The  curves  given  in  Appradix  H  constitute  'ttie  extent  to  which 

we  compared  the  wind  tunnel  model  with  that  arising  out  of  the  real 

data  for  the  coefficient  C^.  Time  permitting,  he  SMLR  technique 

could  have  been  used  to  generate  an  equation  model  for  C^.  just  as 

was  done  for  C  . 

nx 
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6.  SUMMARY  OF  RESULTS  AND  CONCLUSIONS 

The  EBM  methodology  has  been  tested  on  actual  Right  test  data. 

The  Rrst  step  of  the  EBM  method  processed  the  eighteen  maneuvers 

of  data  which  had  been  collected  over  a  period  of  ten  months  and  it 

estimated  an  angle  of  attack  bias  of  2"  and  a  total  airspeed  bias  of 

-3. 9  m/sec  for  the  first  eight  maneuvers  and  none  for  die  ten  maneu- 

vers  of  Right  four.  Scale  factors  of  angle  of  attack  and  sideslip  angle 

were  estimated  as  0.  78  and  0.805.  respectively.  Large  biases  were 

estimated  for  a  number  of  the  measured  variables.  The  roll  and  pitch 

rates  have  the  large  biases  of  -6. 1  deg/sec  and  -2. 7  deg/sec,  respec- 

2 

tively.  The  bias  of  the  fore-aft  acceleration  is  0. 6  m/sec  and  the 

2 

normal  acceleration  (at  c.  g. )  is  -1.  6  m/sec  while  the  normal  acceler- 

2 

ations  off  c.  g.  are  1. 0  m/sec  .  The  bias  of  the  aids  alignment  of  the 
rate  gyro  is  3°.  The  noise  levels  on  the  yaw  rate,  the  rudder  deRection, 
the  total  airspeed  and  the  acceleration  measurements  are  unduly  high 
due  to  the  use  of  a  commutator  to  handle  many  measured  quantities 
rather  than  the  use  of  a  single  channel  for  each  quantity.  There  are 
losses  of  the  airspeed  measurement  over  eirtended  intervals  of  Right. 

The  roll  angle  measurement  was  not  accurate  near  a  roll  angle  of 
180°.  The  vibrations  in  the  nose  boom  caused  large  oscillatory 
deRections  in  the  angle -of- attack  and  sideslip  angle  measurements. 

The  estimation  results  presented  in  Appendix  B  demonstrate  that  the 
first  step  of  the  EBM  method  performed  well  on  this  test  Right  data. 

The  second  step  of  the  EBM  method  identified  detaRed  nonlinear 

models  of  C  .  C  ,  C  and  C  with  partial  results  for  C  and  C  . 

y*  f '  n  m  X  z 

C,  ,  C  and  Crr.  are  identiRed  as  nonlinear 

Ip 

ftmcRons  of  a  and  8  and  the  derivative  Cm  is  identified  as  a  non- 
linear  function  of  a  and  6  .  The  identified  models  and  derivatives  are 
compared  with  the  wind  tunnel  model.  Many  of  the  identified  derivatives 
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match  well  with  the  wind  tunnel  model.  For  example.  Figure  4.2-6 

shows  some  excellent  xnatches  for  with  a  between  18*  aiul  30* 

and  with  |  j3|  between  0*  and  10*.  Fi^ire  5. 2-6  presents  excellent 

matches  for  with  6  between  -28*  and  -24*  and  with  a  between 

m  e 

15*  and  39*. 


The  identified  nonlinear  models  of  the  control  derivatives 
C,  and  C,  are  excellent  znatches  to  die  wind  tunnel  model.  The 

S'  ^ 

a  r 

identified  derivatives  C  and  C  have  a  smaller  magnitude  than 

r  r 

the  wind  tunnel  model  for  some  angles  of  attack. 


The  identified  dynamic  derivatives  are  compared  with  the 
theoretical  prediction  models  of  Ref.  [2].  The  identified  compares 

well  for  all  alpha  but  only  matches  at  low  alpha.  At  hi^  alpha 

C  is  predicted  to  have^values  above  .  10  but  the  identified  model  shows 
P 

that  the  values  are  around  -.075.  The  simulation  study  of  Ref.  [1]  found 

that  the  theoretical  prediction  of  is  inaccurate  at  high  a  and  that 

P 

a  constant  value  of  -.06  provided  a  better  model  for  generating  synthetic 

responses  that  mimics  more  closely  the  real  data.  The  identified  model 

of  obtained  by  processing  the  actual  T-2C  flight  data  agrees  with 
P 

the  modification  of  the  theoretical  model  of  C  obtained  by  trial  and 

n 

error  in  Ref.  [1].  ^ 

The  identified  models  of  and  C ,  agree  fairly  weU  with 

P 

the  theoretical  prediction  models  with  the  exception  that  the  identified 

model  hovers  about  a  zero  value  around  stall  but  the  theoretical  model 

shows  positive  values.  The  identified  model  of  drops  off  more 

r 

sharply  after  stall  than  the  theoretical  does. 
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The  modeling  results  of  Ihis  study  demonstrate  that  die  EBM 
system  identification  method  provides  accurate  nonlinear  modeling  of 
high  a//3  aerodynamics  stability  and  control  characteristics  from 
actual  flight  test  data. 
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7.  RECOMMENDATIONS 

We  recozximend  the  following  fo-'  improving  the  flight  test  data 
for  system  identification  studies: 

1.  Instrumentation  Improvement 

A.  Measure  the  Euler  yaw  angle  0 .  This  will  permit  better 
accuracy  in  estimating  the  bias  on  the  yaw  rate  r. 

B.  Use  sensors  that  measure  the  absolute  values  of  Ihe  Euler 
pitch  and  bank  angles,  6  and  0,  rather  than  measurements 
that  are  relative  to  arbitrary  initial  settings.  This  will  allow 
the  use  of  all  maneuvers  to  estimate  the  biases  of  6  and  0, 
and  it  will  provide  a  check  on  the  condition  d  =  a  at  trim. 

C.  Measure  critical  engine  parameters  for  estimating  thrust 
of  engines.  This  will  permit  using  the  actual  values  of  a  vary¬ 
ing  thrust  to  identify  a  better  drag  model. 

D.  Use  bank  angle  sensors  that  measure  accurately  through 
rolls  greater  than  ISO”.  This  will  permit  better  estimation  re¬ 
sults  at  these  bank  angles  during  high  roll  rates  and  it  will 
provide  continuity  through  0  *  180". 

0  * 

E.  Measure  reliably  the  angular  acceleration  rates  p,  q  and 
r.  This  will  provide  a  check  on  o£f-c.  g.  acceleration  measure¬ 
ments  and  more  accxirate  estimation  of  C .  (t),  C  (t)  and  C  (t). 

jt  m  n 

F.  Use  nose  boom  with  less  vibrational  characteristics  than  is 
on  the  T-2C.  This  will  provide  better  a  and  ^.data  for  estima 
tion  and  modeling. 

G.  Use  airspeed  sensors  which  do  not  give  a  loss  of  measure¬ 
ments  over  extended  intervals  of  flight.  This  permits  better 
estimation  and  modeling  results. 


H.  Use  sin^e  channels  for  all  measured  variables.  This 
provides  the  best  signal- to -noise  ratio  for  a  measured  quantity 
as  con^)sured  to  using  a  commutator. 

I.  Use  pulse  coded  modulation  (PCM)  rather  than  pulse  ampli- 
tude  modulation  (PAM).  This  will  provide  cleaner  signals. 

Rigoiwis  Calibration  and  Consistency  Checkout 

A.  Calibrate  instrumentation  using  static  data  gathered  in 
hangar  (i.  e. ,  with  aircraft  in  hangar,  isolated  from  wind, 
turn  instrumentation  power  on  and  record  a  100  seconds  of 
data).  This  will  permit  the  isolation  and  removal  of  unneces¬ 
sary  biases  in  linear  and  angular  accelerations,  in  p.  q  and 
r.  and  in  6  and  ^ 

B.  Calibrate  instrumentation  using  dynamic  data  gather  on  long 
runway  (i.  e. ,  record  data  with  aircraft  jetting  down  runway  in 
takeoff  phase).  This  will  permit  the  isolation  and  removel  of 
unnecessary  biases  in  a  and 

C.  Calibrate  instrumentation  using  dynamic  data  gather  in 
trim  flight  for  various  trim  values  of  sideslip  angle  and  total 
airspeed.  This  will  permit  the  isolation  and  removal  of  biases 
in  the  controls  for  equilibrium  lli^t. 

D.  Check  out  the  consistency  of  the  measured  variables  using 
several  representative  maneuvers  of  flight  data.  Use  the  equa¬ 
tions  of  motion  to  test  for  consistency.  This  will  provide  a 
checkout  on  scale  factors  and  axis  alignment  as  well  as  the 
compatibility  of  the  measured  variables. 

E.  Fly  calibration  maneuvers  before  and  after  each  flight 
set  of  maneuvers. 


m.  togrove  Design  of  Maneuvers  for  System  Identification 

A.  Fly  more  maneuvers  and  vary  the  variables  such  as 
trim  airspeed  and  sideslip  angle. 

B.  Design  the  maneuvers  so  that  a  good  amount  of  data 
is  provided  at  high  and  low  a  and  at  high  ^  extremes 
rather  than  a  few  isolated  data  points. 

C.  Start  and  end  each  maneuver  with  long  trim  times. 

IV.  Precise  Calculation  of  Aircraft  Constants 

A.  Compute  the  aircraft  mass  for  the  beginning  and  the  end 
of  each  maneuver. 

B.  Compute  the  limits  on  the  changes  in  the  moments  of 
inezHia  due  to  wei^t  distribution  changes  such  as  fuel 
consumption  during  maneuvers. 

C.  Provide  calculations  on  the  location  of  the  center  of 
gravity  with  respect  to  the  aerodynamic  center  for  each 
maneuver. 
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APPENDIX  A 

COMPARISON  OF  AIRSPEED  MEASUREMENTS 

Two  instruznexits  were  used  to  measure  the  total  airspeed  of 
the  T-2C.  One  instrument  had  a  measurement  range  of  500  knots. 

Its  measurement  value  is  called  a  "coarse"  value.  The  other  instru¬ 
ment  had  a  measurement  range  of  100  knots  and  had  the  capability  to 
be  centered  on  any  increment  of  50  knots.  Its  measured  value  is  called 
a  "fine  value". 

The  Figures  A-1  through  A-3  show  a  comparison  between  the 
time  histories  of  the  coarse  and  fine  values.  The  coarse  value  is  denoted 
by  a  solid  line  and  the  fine  value  is  plotted  by  -  dashed  line  with  crosses 
on  it.  In  most  cases  the  dashes  are  so  dose  together  that  the  "fine  value" 
curve  may  appear  as  a  solid  line  with  crosses  on  it.  All  eighteen  maneu¬ 
vers  are  covered  in  these  figures.  Note  that  there  is  an  approximately 
four  meter/ sec  difference  between  the  coarse  and  fine  values  for  all  man¬ 
euvers  except  those  in  flight  number  four  (i.  e. ,  P4M1,  F4M2,  .  .  .  , 
F4M10). 


^=5i 

»p.n" 

a*,t, 


ih 


’•r 


..1f('»»‘">^|u" 

'>'Tii,iiiiii;'!i^i 


iSa-gl 


IlildU 


""'  '  '‘"‘^“^JjiiillilUiifllb:>..i>'d  luiil 


ir*i"T',;.7ii,ii;i 


.•Y' 

*  <  'I  I 

*  ...ll!  .illll 


hi 


.os  lBt.61  t 


iPiffffPl— 1!HI! 

ilia.  :r'‘PHHHi^iiiiMi 


for  maneuvers  FI  Ml  through  F4M1 


iCiailBBSISSBr 

HWiiiipiiHijpiiflijfe: 

. iiisr 


lAi  .  i> 


iilSB 


’  I 


_  i’’"'UL 

iniiii.'. ..... .liiidi 


I ' r'i 


fmm'mffr 


■ii 


a 


I 

I  Fig.  A-2 

i 


ad 

1 


A  comparison  between  the  time  histories  of  the  coarse 
total  airspeed  measurement  for  maneuvers  F4M2  through  F4M9 
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The  estimatioii  results  which  are  an  output  of  the  first  step 
of  the  EBM  method  are  presented  for  all  maneuvers  in  Figures  B-1 
through  B-18.  The  results  cover  four  pages  for  each  maneuver. 

The  first  page  contains  the  angle -of-attack  and  sideslip  results.  The 
second  page  covers  the  time  history  of  the  three  controls:  elevator, 
aileron,  anrf  rudder.  The  third  page  contains  airspeed,  roll  rate, 
pitch  rate,  yaw  rate,  pitch  angle,  roll  angle,  altitude,  and  normal 
acceleration  at  the  center-of-gravity.  The  fourth  page  consists  of 
the  remaining  eight  accelerations:  the  normal  acceleration  at  the 
right  and  left  wing  tips  and  at  the  nose  and  tail,  the  lateral  accelera¬ 
tion  at  center-of-gravity,  nose  and  tail,  and  the  fore  -  aft  acceleration 
at  the  nose. 

These  figures  give  a  comparison  between  the  estimation  results 
of  the  first  step  of  the  EBM  method  and  the  measured  values  of  the 
T-2C  flight  teat  data  that  were  an  input  to  the  EBM  method.  The  esti¬ 
mated  values  satisfy  the  differential  equations  contained  in  Equations 
CA-1)  through  CA-9)  of  Appendix  A  of  Volume  2. 
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APPENDIX  C 


MODELING  RESULTS:  C  MODEL  IN  EQUATION  FORM 

3^ 


This  appendix  contains  the  parametric  model  of  the  side  force 
coefficient  C  ,  identified  using  the  EBM  system  identification  method. 

y 

The  equation  model  identified  for  C  is: 

^  2 

C  »  C  (a)  -6^  +  0  (a)  .  ^  +  L  C  (a) .  f  03.  jS  ^  ) 

y  i'l  Vi 

where  and  are  knots  on  the  l^j  axis  satisfsdng 

and  where  f  03.  ^i.  i'S  a  cubic  spline  function  defined  as 


(■ 


OS|J<S 


f  (-13.  0..  0y^y)  =  -f  03.  0>O 


The  identified  parameters  C  .  C  .  C  .  and  C 


''’o 


are  given  in  Table  C*1  as  functions  of  the  angle  of  attack  a.  The  tabulated 
values  of  the  parameters  hold  at  the  midpoints  of  the  a  intervals  listed  in 
the  first  column  of  the  table.  Interpolation  is  to  be  used  in  between  mid¬ 
points.  The  last  column  in  the  table  lists  the  knot  values  of  the  cubic  spline 
function.  For  example,  at  the  midpoint  a  s  19.5**.  the  knots  are  0^  *  0*. 

3  5a  and  s  18*.  In  order  to  use  the  above  equation  the  function 
2  <3 


must  be  divided  by  the  cube  ot  57. 29  whenever  $,  0^  and 
are  in  degrees. 

In  the  first  column  of  the  table  the  L.  H.  denotes  the  left  hand 
endpoint  and  R.H.  die  right  hand  endpoint. 


Interval 

(Deg) 


R.H 


r 

(per  rad) 


20 

• 

21 

* 

22 

23 

24 

25 

30 

35 

• 

\ 

(per  rad) 


-.661 

-.85 

-.788 

-.908 

-.919 

-.84 


-.781 

-.621 

-.67 

-.76 

-1.08 

-1.03 

-.88 

-.86 

-.922 

-.922 

-.922 

-.922 


Cubic  Spline  Function  Parameters 


knots  on  |^i  axis 
(Deg) 


-2. 

-2. 

-2. 

1.31 

1.54 

-1.33 

-6.98 

-7.80 

-6.58 

-9.2 

I 

-5.8 
1.13  ; 

6. 92  ! 

-5.15 
-4.56  I 
5.49  ^ 


76.56  i 
40.55 
18.54 
-8.58 


0,  5,  18 


0,  5.  20 


0,  5,  24 
0.  5,  25 
0.  5,  28 
0«  5.  28 


C-3 
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APPENDIX  D 


MODELING  RESULTS:  MODEL  IN  EQUATION  FORM 

This  appendix  contains  the  parametric  model  of  the  rolling 

moment  coefficient  C.,  identified  using  the  EBM  system  identification 

X 

method.  The  equation  model  identified  for  is: 


(a)*6^  +  C^  (a) -6^  +  0^  (a).j3+  L  (a) •  f 03.  )3.^j) 


-‘p  ^ " 


C, 

2V 


where  and  are  knots  on  the  l0 1  axis  satisfying  0  s  g  <  ^2  ^  ^3  ^  ^4 
and  where  f  (fl,  is  a  cubic  spline  function  defined  as 


f(^.  ^i.  ^i^i)  >  \  i 


0^  fi. 


S  ^6^  B 
*1  ^  '^i+l 


f  (-P.  * 


-f  O.  P>0 


The  identified  parameters  .  i  *  1,  2.  3 . 

6  6/3^ 

a  r  o  /T 

and  are  given  as  functions  of  the  angle  of  attack  in  Tables  D-1 

P  ** 

and  D-2.  The  tabulated  values  of  the  parameters  hold  at  the  midpoints  of 
the  a  intervals  listed  in  the  first  coluznns  of  these  two  tables.  The  knots 
of  the  cubic  spline  function  are  given  in  the  last  column  of  Table  D-1.  The 
function  f  O,  must  be  divided  by  the  cube  of  57. 29  whenever 

fi,  and  are  in  degrees. 


D-1 


Table  D-1 


Identified  parameters  of 


Cubic  Spline  Function  Parameters 


C. 


.011 

.0076 

.0041 

.0113 

.0097 

.0065 

.0072 

.0066 

.0061 

.0135 

.0127 

.012 

.0145 

.0190 

.0093 

.0123 

.0144 

.0012 


UW'M 


.0068 

.0068 

.0068 

.0068 

.0068 


(per  rad  ) 


-.484 

-.484 

1.33 

1.05 

1.17 

.78 

.43 

.21 

-1.141 

-3.05 

14.23 

20.41 

-9.44 

-2.01 

-.66 


53  -10.83 


-6.9 


1.33 

.39 

-1.43 

-1.38 

-1.58 


" 

1 

i 

0, 

0. 

0. 

10 

10 

12 

1 

1 

0. 

13 

1 

0, 

13 

- 

0, 

14 

- 

0, 

14 

- 

0. 

14 

-3.44 

0. 

5,  10, 

16 

-4.79 

0. 

5,  10, 

16 

102.2 

0, 

2,  10, 

15 

97.17 

0. 

2,  10, 

15 

0. 

10,  15 

- 

0, 

10,  18 

■ 

0, 

10,  18 

0. 

5,  18 

0. 

5,  18 

- 

0, 

5,  20 

- 

0. 

5,  20 

0, 

5,  23 

- 

0. 

5,  24 

- 

0. 

5,  25 

0, 

5,  28 

- 

0, 

5,  28 

a 

Interval 

(Deg) 


Table  D-2 

Identified  dynamic  derivatives  of 


P 

(per  rad) 


-.6100 

-.6100 

-.5957 

-.5861 

-.5766 

-.5670 

-.5460 

-.5440 

-.5010 

-.4940 

-.4400 

-.4680 


A  I 

18 

lO  1  • 

19  i 

19 

20  j 

20 

21 

21 

22  j 

22 

23 

23 

24 

24 

25 

25 

30 

30 

35 

35 

40  1 

(per  rad) 


.1102 

.1467 

.1613 

.1759 

.1904 

.205 

.2196 

.2342 

.2488 

.2634 

.278 

.290 

.288 

.312 

.326 

.310 

.357 

.183 

.218 

.232 

.162 

.116 

.129 

.1261 

.1232 

.1202 

.1115 

.0969 

.0823 


APPENDIX  E 


MODELING  RESULTS:  MODEL  IN  EQUATION  FORM 


This  appendix  contains  the  parametric  model  of  the  yawing 

moment  coefficient  identified  using  the  EBM  system  identification 

method.  The  equation  model  identified  for  C  is: 

n 

3 

C  »C  (a)-6  +C  (a)-/5+  E  C  (a) •  f (ji  3., 
n  n  r  n«  n  ,  i  i+i 

*r  Po 

<=n 

P  r 

where  and  are  knots  on  the  |^|  axis  satisfying  0  ^ 

and  where  f  (fi,  0^,  0^^^)  is  a  cubic  spline  function  defined  as 


f  (0,  Hj, 


/  0 


j 


A  ^0^  0 
^  ^  '^i+l 


^i+1  ~  "  3  ^i  "  3  ^i+1^ 


0  <0 
^i+1  ^ 


U-0,  0^,  =  -i{0,  0^^^)  0>O 


The  identified  parameters  C  .  C  ,  C  «  i  -  1,  2,  3, 

“6  “/5  “«3 

r  o 

C  and  C  are  given  as  functions  of  the  angle  of  attack  in  Tables  E-1 
n  n 

P  r 

and  E-2.  The  tabulated  values  of  the  parameters  hold  at  the  midpoints 
of  the  a  intervals  listed  in  the  first  columns  of  these  two  tables.  The 
knots  of  the  cubic  spline  function  are  presented  in  the  last  column  of 
Table  E-1.  The  function  f(A  ^^+1^  must  be  divided  by  the  cube  of 


57.29  whenever  and  are  in  degrees. 

In  the  first  column  of  Tables  E-1  and  E-2  the  L.H.  denotes 
the  left  hand  endpoint  of  the  cl  interval  and  R.H.  the  right  hand  exidpoint. 


I 


Table  E-1 


ldentifi.ed  parameters  of  C 


a 

Interval 

(Deff) 

C 

r 

rr 
!  % 

(per  rad) 

C 

“  3 

d 

^  3 

(per  rad  ) 

(per  rad 

L.H. 

R.H. 

(per  rad) 

-4 

2 

-.044 

.076 

.53 

- 

2 

3 

-.044 

.065 

1.03 

- 

3 

4 

-.044 

.051 

1.63 

- 

4 

5 

-.044 

.051 

1.63 

- 

5 

6 

-.044 

.064 

1.63 

- 

6 

7 

-.044 

.073 

.55 

- 

7 

8 

-.044 

.096 

-.473 

- 

8 

9 

-.041 

.089 

-.397 

- 

9 

10 

-.038 

.081 

.02 

- 

10 

11 

-.040 

.091 

-1.03 

- 

11 

12 

-.031 

.082 

-.42 

- 

12 

13 

-.027 

.058 

4.93 

-17.33 

13 

14 

-.032 

.052 

4.18 

-15.93 

14 

15 

-.030 

-.013 

89.74 

-6.91 

15 

16 

-.022 

.032 

-5.35 

.54 

16 

17 

-.017 

.011 

.76 

-15.4 

17 

18 

-.025 

.015 

-1.86 

9.45 

18 

19 

-.020 

-.023 

3.15 

-6.32 

19 

20 

-.014 

-.016 

.43 

1.73 

20 

21 

-.011 

-.011 

.79 

-.29 

21 

22 

-.014 

.005 

-5.58 

2.46 

22 

23 

-.015 

-.013 

-4.32 

1.84 

23 

24 

-.014 

-.015 

.29 

.29 

24 

25 

-.016 

.006 

• 

o 

1-4 

1 

o 

25 

30 

-.017 

-.011 

.46 

1 

.46 

30 

35 

-.019 

-.020 

.55 

.55 

35 

40 

-.021 

.034 

.80 

.80 

“  3 

B 

^  3 

(per  rad  ) 

Knots  on  \B  \ 
axis 
(Deg) 

- 

0.  10 

- 

0,  10 

j 

0«  10 

- 

0,  10 

1 

0,  10 

- 

0,  12 

- 

0.  13 

- 

0,  13 

- 

0,  14 

- 

0,  14 

- 

0.  14 

38.52 

0.  16 

32.19 

0.  5.  10,  16 

-1.55 

0,  5,  10,  16 

-20.38 

0,  2,  10,  15 

- 

0,  2.  10,  15 

- 

0,  10,  15 

- 

0,  10,  18 

- 

0,  5,  18 

- 

0,  5,  18 

-  1 

0,  5,  20 

1 

0,  5,  20 

^  i 

0,  5.  23 

- 

0,  5,  24 

i 

0,  5,  25 

- 

0,  5,  28 

- 

0,  5,  28 

APPENDIX  F 


MODELING  RESULTS:  C  CURVES 

m 


A  comparison  between  the  wind  tunnel  values  of  and  the  values  of 

extracted  from  the  actual  T-2C  flight  test  data  is  contained  in  Figures 

F-1  through  F-14.  In  these  figures  the  aerodynamic  coefficient  C  is 

m 

shown  as  a  function  of  the  angle  of  attack  for  fixed  one  "ee  intervals 
of  the  elevator  control.  The  range  of  the  elevator  control  6^  covered 
is  from  -28°  to  0°. 

Figures  F-1  through  F-14  were  generated  as  follows:  a  particularly 

o 

small  range  of  the  elevator  control  was  selected  from  8  to  6  +1  where 

6  G 

8  is  an  integer  between  -28°  and  -1°.  All  C  (t)  estimates  (i.  e. ,  an 
e  m 

output  of  the  first  step  of  the  EBM  method)  were  collected  together  from 

all  maneuvers.  From  this  collection  those  C  (t)  estimates  having  a 

m 

corresponding  elevator  control  8  (t)  within  the  selected  range  8  to  8  +1° 

e  e  e 

were  ordered  according  to  the  angle  of  attack  variable  and  then  plotted 

as  one  of  the  "real  data  values"  curves  in  Figwe  F-1  through  F-14.  The 

"wind  tuimel  values"  curves  were  obtained  by  evaluating  the  wind  tunnel 

model  using  the  estimated  states  and  the  elevator  control  8^(t).  This  gave 

a  collection  of  wind  tunnel  coefficient  values  WTC  (t)  from  all  maneuvers. 

m 

From  this  collection  those  WTC^(t)  values  having  a  corresponding 

elevator  control  8  (t)  within  the  selected  range  8  to  8  +1  were  ordered 
e  ®  e  e 

according  to  the  angle  of  attack  variable  and  then  plotted  as  one  of  the 

"wind  tunnel  values"  curves  in  Figure  F-1  through  F-14.  It  should  be 

observed  that  the  curves  in  these  figures  have  a  dependency  on  $,  q 

and  if  built  into  them  by  the  nature  of  the  maneuvers. 
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MODELING  RESULTS;  C  CURVES 
_ z _ 

A  comparison  between  the  wind  tunnel  values  of  C  and  the 

z 

values  of  extracted  from  the  actual  T'2C  flight  test  data  is  contained 

in  Figures  G-1  through  G-14.  In  these  figures  the  aerodynamic  coefficient 

C  is  shown  as  a  function  of  the  angle -of -attack  for  fixed  one  degree  in- 
z 

tervals  of  the  elevator  control.  The  range  of  the  elevator  control  covered 
is  from  -28°  to  0°. 

Figures  G-1  through  G-14  were  generated  as  follows:  a  particularly 

small  range  of  the  elevator  control  was  selected  from  6  to  6  +  1° 

e  e 

where  6  is  an  integer  between -28  and -1.  All  C  (t)  estimates  (i.e., 

€  Z 

an  output  of  the  first  step  of  the  EBM  method)  fiom  all  maneuvers  haviiig 

a  corresponding  elevator  control  6  (t)  within  the  selected  range  6  to 

o  ^  ^ 

6^  +  1  were  ordered  according  to  the  angle -of -attack  variable  and  then 

plotted  as  one  of  the  "real  data  values”  curves  in  Figures  G-1  through 

G-14.  The  "wind  tunnel  values"  curves  were  obtained  by  evaluating  the 

wind  tunnel  model  using  the  estimated  states  which  have  an  associated 

elevator  control  6  (t)  in  the  selected  range  6  to  6  +  1°,  by  storing 

the  resulting  wind  tunnel  coefficient  values  WT  C  (t),  by  ordering  them 

z 

with  respect  to  the  angle -of -attack  variable.*,  and  by  plotting  the  ordered 
wind  tuimel  values.  It  should  be  observed  that  the  curves  in  Figures  G-1 
through  G-14  have  a  dependency  on  )3,  q  and  a  built  into  them  by  the 
nature  of  the  maneuvers. 
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MODELING  RESULTS:  C  CURVES 

X 


A  comparison  between  the  wind  tunnel  values  of  C  and  tiie 
values  of  extracted  from  the  actual  T-2C  flight  test  data  is  contained 
in  Figures  H-1  through  H-14.  The  coefficient  is  given  as  a  function 
of  the  angle-of-attack  for  fixed  one  degree  intervals  of  the  elevator  control. 
The  range  of  the  elevator  control  covered  is  from  -28®  to  0®. 


The  Figures  H-1  through  H-14  were  obtained  in  a  way  similar  to 
those  contained  in  Appendix  G.  A  discussion  is  given  therein. 
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APPENDIX  I 


MODELING  RESULTS:  C  MODEL  IN  EQUATION  FORM 

TO. 


This  appendix  contains  the  parametric  model  of  the  pitching 

moment  coefficient  C  ,  identified  using  the  EBM  system  identification 

method.  The  equation  model  identified  for  C  is: 

m 

6 

•  (a  -  a  )  +  E  C  (6  )  •  g  (a,  a  ,  a  ) 
m  me  m_  e  i  m »  e  i  i+i 

O  1*1  jL 

o 


+  C  (a)  .  q  +  <a)  .  a  +  E  (Q  •  h  (/5. 

“q  “d  3*1  “^2  J 


where  a . . . ,  are  the  knots  on  the  a-axis  satisfying  -4  <  < 

—  <  a^,  where  fiy  are  the  knots  on  the  1^  ]  -  axis  satis^^ing 

0  s  /9^  <  where  g  (a,  a^,  is  a  quadratic  spline  function 

defined  as 


g  (a. 


V  ®i+l^ 


0  a<a. 

1 


Vj  -  “i*  -  T“i  -  i  Vl> 


Vi^“ 


and  where  h  {fi,  fiy 


is  a  quadratic  spline  function  defined  as 


0  1)31  <  ^  j 


h  V.  fy 


The  identified  parameters  C 


used  in  between  midpoints.  The  knots  are  presented 

in  the  table  for  each  6  interval.  For  example,  at  6  =-14.5%  the 

G  0 

knots  are  =  14®,  *  15®,  =  18®,  Og  =  20®  and  =  24®. 

The  function  g  (a,  a^,  must  be  divided  by  the  square  of  57. 29 

whenever  a,  a^,  are  in  degrees. 


The  identified  parameters  C  ,  i  =  1.  2,  3.  C  ,  C  are 

02  q  a 

given  in  Tables  1-2  and  1-3  as  functions^of  the  angle  of  attack  a.  The 
tabulated  values  of  these  parameters  hold  at  the  midpoints  of  the  a  in¬ 
terval  listed  in  the  first  column  of  the  table.  Interpolation  is  to  be  used 
in  between  midpoints.  The  knots  0^,  0^,  0^,  0^  are  given  in  Table  1-2 
for  each  a  interval.  For  example,  at  a=  14.5®,  the  knots  are  0^  =  0®, 
02  ~  0^  -  0^  *  15**  The  function  h  {0  0^,  must  be  divided 

by  the  square  of  57.29  whenever  are  in  degrees.  ' 


Table  I-l 

Identified  a>parameters  of  C 

xn 


6 

Interval 

(Deg) 

L.H.  R.H. 

o 

^m 

(per  rad) 

a  Range 
_ (Deg)  _ 

Low 

High 

-28 

-27 

.0993 

-4.51 

34.82 

-27 

-26 

.1152 

-4.78 

37.54 

-26 

-25 

.084 

-4.71 

15.56 

37.51 

-25 

-24 

.066 

-3.59 

16.09 

36.99 

-24 

-23 

.067 

-1.90 

14.00 

32.76 

-23 

-22 

.060 

-2.13 

14.00 

31.12 

-22 

-21 

.007 

-.51 

14.00 

29.97 

-21 

-20 

.117 

-5.13 

14.00 

28.77 

-20 

-19 

.256 

-1.54 

7.17 

27.55 

-19 

-18 

.219 

-1.03 

7.80 

26.61 

-18 

-17 

.204 

-1.18 

8.08 

25.34 

-17 

-16 

.145 

.66 

25.38 

-16 

-15 

.168 

.04 

4.04 

24.66 

-15 

-14 

.164 

.68 

1.75 

23.69 

-14 

-13 

.157 

.33 

1.65 

23.02 

-13 

-12 

.158 

.183 

1.42 

23.40 

-12 

-11 

.156 

-.117 

1.01 

23.49 

-11 

-10 

.144 

.008 

23.43 

-10 

-9 

.112 

.026 

22.51 

-9 

-8 

.107 

-.115 

22.01 

-8 

-7 

.133 

-.415 

-1.66 

20.57 

-7 

-6 

.150 

-.501 

-3.00 

20.27 

-6 

-5 

.109 

-.21 

-3.00 

19.70 

-5 

-4 

.143 

-.83 

-3.00 

19.99 

-4 

-3 

.143 

-.98 

-3.00 

20.00 

-3 

-2 

.125 

-.795 

-3.00 

18.83 

-2 

1  -1 

.111 

-.92 

-2.00 

18.93 

-1 

mm 

.072 

-1.03 

-1.00 

18.94 

Model 

Quadratic  Spline  Function 

Parameters _ 

Knots  on  a  axis 
_ (Deg) _ 

15«  18,  25,  30,  35 

15,  18,  25.  30,  39 

15,  18.  25,  30.  39 

15,  18.  25,  30,  39 

14.  15.  18.  20,  25,  35 

14.  15.  18.  20.  25.  33 

14.  15,  18.  20,  25.  32 

14,  15.  18,  20.  25.  31 

7.  14.  15.  18.  20.  25.  29 

7.  14,  15.  18,  20,  25,  28 

7.  14,  15,  18,  20,  25.  27 

5,  14.  15,  18.  20.  25.  26 

4.  14,  15.  18,  20,  25 

1,  14,  15,  18.  20.  24 

1.  14,  15,  18.  20.  24 

1.  14.  15.  18.  20.  25 

1,  14.  15.  18.  20,  25 

1.  14.  15.  18.  20.  25 

1.  14.  15.  18.  20.  24 

1.  14.  15.  18.  20.  23 

-2,  14.  15,  18,  20.  22 

-3,  14.  15,  18.  20,  21 

-3.  14.  IS.  18.  20 

-3,  14,  15,  18,  20 

-3,  14.  15,  18,  20 

-3,  14.  15,  18,  19 

-2.  14.  15.  18.  19 

•1.  14.  15.  18.  19 


Table  I-l  (cont'd) 


1  Quadratic  Spline  Function  Parameters  (per  rad^)  | 

Interval 

(Deg) 

L.H.  R.H. 

2 

C‘ 

2 

c 

“  2 
a 

4 

C 

C 

“  2 
®6 

-28 

-27 

85.76 

-7.34 

11.21 

-3.09 

- 

- 

-27 

-26 

85.76 

-7.34 

11.21 

-3.09 

- 

- 

-26 

-25 

85.76 

-7.34 

11.21 

-3.09 

- 

- 

-25 

-24 

60.55 

-7.68 

14.72 

-3.12 

- 

- 

-24 

-23 

0.0 

13.90 

9.63 

15.87 

-19.15 

- 

-23 

-22 

0.0 

28.85 

-32.18 

29.14 

-21.8 

- 

-22 

-21 

0.0 

-17.37 

3.23 

21.44 

-21.32 

- 

-21 

-20 

69.02 

0.0 

-27.3 

28.14 

11.57 

-17.58 

-20 

-19 

-8.45 

103.88 

-13.05 

17.48 

12.49 

-21.17 

-19 

-18 

-12.66 

82.4 

1.96 

12.21 

6.37 

-18.92 

-18 

-17 

-9.09 

77.8 

-11.54 

36.96 

-2.63 

17.72 

-17 

-16 

-20.87 

89.72 

-14.22 

58.31 

-16.66 

26.54 

-16 

-15 

-12.37 

23.17 

9.69 

29.43 

-11.48 

- 

-15 

-14 

-13.6 

54.2 

1.30 

63.41 

-55.5 

- 

-14 

-13 

-10.83 

-8.55 

38.94 

27.3 

-68.6 

- 

-13 

-12 

-10.60 

30.2 

27.84 

22.4 

-50.8 

- 

-12 

-11 

-8.U 

-13.65 

55.05 

-40.84 

-6.91 

- 

-11 

-10 

-9.56 

73.9 

-5.70 

48.6 

-29.8 

- 

-10 

-9 

-8.48 

21.83 

21.84 

10.37 

-27.75 

- 

-9 

-8 

-7.28 

-45.99 

60.42 

-45.77 

34.4 

- 

-8 

-7 

-2.95 

-17.85 

2.94 

56.71 

-13.75 

- 

-7 

-6 

-2.61 

-20.18 

9.83 

43.06 

- 

- 

-6 

-5 

-4.28 

9.96 

11.87 

-70.1 

- 

- 

-5 

-4 

-.69 

-152.6 

89.53 

-204.1 

- 

- 

-4 

-3 

-.11 

■m.i 

-16.75 

148.0 

- 

- 

-3 

-2 

-1.62 

-41.6 

23.57 

176. 

- 

- 

-2 

-1 

-1.24 

-75.0 

28.9 

182. 

- 

- 

-1 

0 

-.95 

-42.5 

32.0 

32.0 

- 

- 

0 

1 

-6.50 

-126.15 

27.58 

- 

- 

1-4 


Table  1-2 

Identified  /3 -parameters  of  C 


(per  rad) 


-.183 


Quadratic  Spline  Function  Parameter  (per  rad^) 


Knots  on  |/9|  axis 
(Deg) 


-.413 
-3.09  ‘ 
-3.14 
-2.68 
-1.06 
-.366 
-10.01 
-10.90 
-12.91 
-6.66 


1.36 


3.18 


13.51 


-3.97 

3.57 

-3.89 


-19.6 


-17.34 

-8.09 

11.47 


0,  15.  10.  16 
0.  5.  10.  16 
0.  5.  10.  15 
0.  5.  10.  15 
0.  5.  10.  15 
0.  10.  18 
0.  10.  18 
0.  10.  18 


0,  10.  20 
0.  10.  20 
0.  10,  23 


Table  1-3 


Identified  dynamic  derivatives  of  C  model 


0 

1 

1 

2 

2 

3 

3 

4 

4 

5 

5 

6 

6 

7 

7 

8 

8 

9 

9 

10 

10 

11 

11 

12 

12 

13 

13 

14 

14 

15 

15 

16 

16 

17 

17 

18 

18 

19 

19 

20 

20 

21 

21 

22 

22 

23 

23 

24 

24 

25 

25 

30 

30 

35 

36 

40 

